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AY. 


ABSTRACT 


A  light  scattering  apparatus  to  measure  particle  size 
( D ^ 2 ^  i-n  a  solid  rocket  motor  was  improved.  Multiple 
consecutive  scans  of  two  photodiode  arrays  were  accomplished 
with  a  pacing  circuit  and  added  memory.  The  device  was 
calibrated  using  various  suspended  particle  samples  and 
found  to  make  accurate  measurements. 
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INTRODUCTION 


I  . 

Performance  prediction  codes  for  solid  rocket  motors 
model  two  phase  flow  losses  as  functions  of  particle  size. 

In  addition,  particle  size  within  the  grain  port  strongly 
affects  the  damping  of  combustion  pressure  oscillations.  At 
present  these  models  are  based  on  particle  size  data  from 
collected  exhaust  samples  [Ref.  1].  However,  particle  size 
varies  with  position  in  the  motor  and  otr.er  parameters 
(pressure,  propellant  formulation,  nozzle  design,  etc.). 
Therefore,  experiments  to  determine  how  particle  size  varies 
in  the  actual  flow  environment  of  the  motor  (i.e.,  across  the 
nozzle)  are  needed  to  validate  the  models  for  two  phase  flow 
losses.  Cramer  [Ref.  2]  and  Karagounis  [Ref.  3]  provide  a 
good  summary  of  the  subject  and  the  Naval  Postgraduate  School 
Combustion  Laboratory  effort  to  obtain  particle  size  data 
across  the  exhaust  nozzle  of  a  solid  rocket  motor. 

The  method  used  in  this  continuing  effort  was  the  dif- 
fractively  scattered  light  technique.  The  diffraction  pattern 
of  light  scattered  by  particles  are  analyzed  to  determine  the 
volume  to  surface  mean  diameter  [Refs.  4  through  11].  This 
method  has  the  disadvantages  that  size  distributions  cannot 
be  easily  determined  and  particles  larger  than  some  threshold 
size  will  not  be  detected  due  to  the  exceedingly  small  angles 
as  which  they  scatter  light.  However,  it  has  the  advantage 


that  it  is  non- intrusive 


ar.c  ir.  z-.e'rv  can  be  used  in  tue 
internal  motor  environment. 

Use  of  this  method  for  particle  sizing  in  solid  rocket 
motors  at  the  Naval  Postgraduate  School  Combustion  Laboratory 
was  begun  by  Karagounis  [Ref.  3].  The  apparatus  was  sub¬ 
sequently  redesigned  and  the  data  acquisition  equipment  up¬ 
graded  with  the  introduction  of  the  Hewlett  Packard  3054A 
data  acquisition  system  with  an  HP  85  as  the  controlling 
computer  [Ref.  12].  The  investigation  by  Cramer  and  Hansen 


followed  and  showed  that  propellent  composition  can  limit 
the  application  of  the  technique.  Large  particulate  combustion 
products  in  the  flow  made  particle  size  data  difficult  to 
obtain.  This  was  especially  true  if  only  one  measurement  of 
the  scattering  profile  was  made  during  a  test  firing. 

To  address  this  problem  in  the  present  study  several 
improvements  were  made.  A  cleaner  burning  propellant  was 
obtained  to  reduce  char  agglomerates  in  the  exhaust  products. 

A  more  statistically  valid  data  sample  (multiple  measurements 
during  a  single  test)  was  made  possible  with  added  momory 
in  the  data  acquisition  equipment  and  a  pacing  circuit  which 
allowed  full  use  of  this  memory.  Data  reduction  was  also 
improved  with  a  Hewlett  Packard  9836S  computer  combined  with 
a  more  recently  developed  approach  to  particle  sizing 
presented  by  Buchele  [Ref.  13].  This  method  is  discussed 


later  in  depth. 


The  focus  of  this  cr.esis  pro] ecu  was  the  following; 

(1)  Implement  the  HP  9836S  as  the  system  controller. 

(2)  Expand  the  multiprcgrammer  memory  in  order  to  obtain  up 
to  eight  consecutive  scans  of  the  diode  arrays  during 

a  test  firing. 

(3)  Improve  data  reduction  techniques  by  the  method  of 
Buchele  [Ref.  13]. 

(4)  Modify  the  apparatus  and  experimental  procedures  to 
improve  the  angular  resolution  and  to  reduce  extraneous 
light . 


(5)  Certify  the  proper  functioning  and  accuracy  of  the 
apparatus  prior  to  actual  motor  testing. 


II.  THEORETICAL  BACKGROUND 


A.  GENERAL  DISCUSSION 

The  completely  general  theory  of  scattering  was  developed 
by  Mie  and  is  presented  by  Van  de  Hulst  [Ref.  14] .  The 
light  scattering  characteristics  for  spherical  particles  of 
any  size  are  fully  described  in  a  mathematical  format.  The 
Mie  scattering  functions  contain  Legendre  polynomials  and 
spherical  Bessel  functions  and  fully  treat  the  phenomena  of 
reflection,  refraction,  diffraction,  and  extinction.  The 
full  theory  is  most  often  applied  when  particle  size  is 
approximately  the  same  as  the  wavelength  of  the  incident  light. 
Van  de  Hulst  [Ref.  14]  calls  this  the  regime  of  Anomalous 
Scattering . 

For  particle  sizes  much  smaller  than  the  wavelength  of 
light  the  Mie  equations  simplify  to  a  form  which  is  more 
dependent  on  the  index  of  refraction  of  the  particles  and 
less  dependent  on  particle  size.  This  is  called  Rayleigh 
Scattering . 

The  study  of  particle  size  behavior  in  solid  propellant 
rockets  mainly  covers  sizes  much  greater  than  the  wavelength 
of  light.  Scattering  by  large  particles  such  as  these  is 
described  adequately  by  Fraunhofer  diffraction. 
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3.  APPLICATION7  TO  LARGER  PARTICLES 

The  ringed  diffraction  pattern  generated  by  a  hole  in 
a  mask,  or  a  number  of  particles  of  the  same  size  is  de¬ 
scribed  by  the  equation: 

2J1 (a6)  2 

1(9)  =  [ — ± - ] 

010 

where: 

1(6)  describes  the  relative  intensity  of  the  scattered 

light  at  an  angle  theta  (8) 

J^(a0)  is  the  Bessel  function  of  the  first  kind 

a  =  is  the  particle  size  parameter  for  diameter 

D  and  wavelength  of  light  lambda  (X) . 

Measuring  the  particle  size  for  a  monodispersion  can  be 
accomplished  by  measuring  the  angular  position  of  a  dark  or 
bright  ring  in  the  diffraction  pattern.  For  a  dark  ring  the 
zero  of  the  Bessel  function  corresponding  to  the  ring  is  set 
equal  to  (a0)  and  particle  size  is  determined  directly  [Ref.  5] 
For  bright  rings  one  sets  (a0)  equal  to  the  corresponding 
maximum  of  the  Bessel  function  and  solves  for  the  diameter. 

The  above  method  is  not  used  for  polydispersions  since 
the  discrete  rings  are  not  observed.  However,  Dobbins,  et  al . 
[Ref.  5]  introduced  a  significant  improvement  in  the  dif- 
fractively  scattered  light  method  of  particle  sizing.  They 
found  that  although  the  method  was  not  directly  able  to 


determine  distributions  of  sizes,  the  volume  to  surface 
mean  diameter  defined  by 


/Q«Nr(D)D3dD 

/Q®Nr(D)D2dD 


(1) 


where: 

Nr(D)  is  a  distribution  function  describing  the  pro¬ 
portion  of  particles  with  diameter  (D)  in  the  sample, 
could  be  accurately  measured, 

A  curve  for  sizing  polydispersions  was  presented  which  was 
used  by  Cramer  and  Hansen  [Refs.  2,  12] . 

Two  phase  flow  losses  are  often  calculated  in  terms  of 
D43>  If  the  distribution  of  sizes  in  the  polydispersion  is 
well  behaved  then  D^2  and  other  diameters  such  as  D43  can 
be  easily  related  [Ref.  15].  Reference  5  reported  that  very 
small  particles  in  the  distribution  have  a  minor  influence 
on  the  scattering  profile.  This  makes  the  measurement 
technique  promising  for  the  two  phase  flow  loss  study 
since  very  small  particles  do  not  contribute  significantly 
to  these  losses  and  so  are  of  less  interest. 

Roberts  and  Webb  [Ref.  6]  essentially  confirmed  the 
conclusions  of  Reference  5  and  presented  a  similar  curve  for 
use  in  sizing. 

More  recently,  Buchele  [Ref.  13]  gives  a  good  summary  of 


experimental  techniques  for  particle  sizing  by  measuring 


diffracted  light.  One  point,  of  interest  in  his  report  is 
that  he  represents  the  scattering  profile  of  a  polydispersion 
with  a  function  which  closely  approximates  the  curves  of 
References  5  and  6. 

I  (6)  =  EXP  - { . 57a6 ) 2 


This  function  from  Reference  13  and  the  curve  from 
Reference  5  were  both  used  in  the  present  study  to  evaluate 
the  apparatus  to  be  used  with  solid  propellant  rocket  motors. 

An  additional  detail  of  measuring  scattering  profiles  is 
covered  by  Van  de  Hulst  [Ref.  14].  The  wavelength  of  light 
used  in  the  scattering  calculations  depends  on  the  index  of 
refraction  of  the  medium  containing  the  particles.  The 
wavelength  used  in  all  calculations  must  be: 


A 


where : 

Xq  is  the  wavelength  of  light  in  a  vacuum  and, 

M  is  the  refractive  index  of  the  medium  with  respect 
to  a  vacuum. 

Thus,  the  size  parameter  (a)  becomes: 


•RDM 


and  the  beam  spread  parameter  becomes: 


Another  consideration  is  presented  by  Gumprecht  and 
Sliepcevich  [Ref.  4].  Light  scattered  by  particles  in  a 
medium  is  refracted  as  it  crosses  each  interface  of  the 
container  holding  the  medium.  This  is  discussed  further  in 
the  section  on  calibration  and  evaluation  of  the  apparatus. 

Additional  complications  arise  with  the  full  treatment 
of  the  index  of  refraction  of  the  particles  with  respect  to 
the  medium.  But,  for  Fraunhofer  diffraction  alone  this  aspect 
can  be  neglected. 

C.  RESTRICTIONS  AND  SOURCES  OF  ERROR 

Some  restrictions  on  use  of  the  method  are  described  by 
Dobbins,  et  al.  [Ref.  5]  and  were  satisfied  as  described  by 
Cramer  [Ref.  2].  These  are  related  to  the  size  of  particles, 
the  distance  to  the  detector,  and  some  phenomenon  covered  in 
the  rigorous  Mie  theory. 

One  must  keep  in  mind  also  that  the  curves  developed  for 
polydispersions  are  based  on  the  Upper  Limit  Distribution 
Function  of  Mugele  and  Evans  [Ref.  7].  This  means  that  no 
particles  with  size  greater  than  approximately  ten  times  the 
mean  should  be  in  the  sample  [Ref.  13].  This  appears  to  be 
a  mild  restriction.  Van  de  Hulst  [Ref.  14]  describes  the 
criteria  for  single  scattering  and  a  simple  test  to  verify 


proportional  to  the  number  of  particles  the  mathematics 
remain  simple. 

Sources  of  error  of  the  dif fractively  scattered  light 
method  are  covered  by  Buchele  fRef.  13]  and  are  presented 
here . 

(1)  Inaccuracy  of  angular  measurement  or  the  limited  ability 
to  resolve  small  angles  and, 

(2)  Inaccuracy  of  the  intensity  measurement  due  to  extraneous 
light . 

Extraneous  light  includes  all  light  other  than  scattered 
light  from  the  particles.  Some  examples  are  scattering  from 
an  aperture  or  dirty  test  section  windows.  Refraction  of  the 
beam  due  to  gas  density  gradients  and  image  point  broadening 
from  turbulence  are  others.  Laser  speckle  is  also  extraneous 
light. 

The  sources  of  error  addressed  in  this  investigation  are 
discussed  in  the  related  portions  of  the  paper. 


III. 


EXPERIMENTAL  APPARATUS 


A  photograph  of  the  apparatus  is  presented  in  Figure  1. 

A  schematic  is  presented  in  Figure  2.  The  light  scattering 
equipment  was  mounted  on  two  optical  benches.  Components  for 
measurements  in  the  exhaust  plane  were  mounted  on  one  bench. 
The  other  bench  held  the  equipment  associated  with  the  motor 
cavity.  The  light  source  was  an  eight  (8)  milliwatt  Helium 
Neon  laser  mounted  on  the  exhaust  bench.  A  collimated  beam 
was  required  so  a  spatial  filter/collimator  was  used.  A 
modification  to  this  collimator  is  discussed  later  in  this 
section.  The  collimated  beam  passed  through  a  cube  beam 
splitter  and  the  second  beam  was  diverted  to  a  90  degree 
prism  on  the  other  bench.  The  original  beam  continued  through 
the  motor  exhaust  plane.  The  other  beam  was  routed  through 
the  nitrogen-purged  glass  windows  in  the  motor  housing. 

Each  beam  was  then  intercepted  by  a  physical  stop  located 
in  front  of  its  set  of  receiving  optics.  The  further  the 
stop  was  placed  from  the  test  section,  the  smaller  the  angle 
at  which  scattered  light  could  be  measured.  In  this 
apparatus,  the  stops  were  placed  approximately  30.5 
centimeters  from  the  exit  plane  of  the  test  section.  This 
allowed  a  minimum  angle  of  approximately  .008  radians  to  be 
measured.  Light  scattered  at  angles  greater  than  this  was 


not 


intercepted  and  continued  past  the  edges  of  the  stcc. 

The  stop  served  to  keep  the  transmitted  beam  out  of  the 
measuring  optics  and  thus  reduce  extraneous  light.  The 
stops  also  improved  optics  alignment.  This  is  discussed 
under  calibration  and  evaluation  of  the  apparatus. 

The  scattered  light  passed  through  a  narrow  pass  filter 
which  admitted  only  light  of  the  Helium  Neon  frequency. 

This  filter  served  to  reduce  extraneous  light  from  the  external 
surroundings . 

An  objective  lens  of  50  centimeter  focal  length  was 
located  behind  the  narrow  pass  filter.  This  lens  imaged 
onto  a  photodiode  array  the  scattering  profile  of  the  particles 
in  the  test  section.  The  shadow  of  the  beam  stop  was  also 
imaged  since  the  stop  was  between  the  test  section  and  the 
objective  lens.  This  was  a  limitation  which  is  discussed 
under  calibration  and  evaluation. 

The  photodiode  arrays  were  the  same  units  used  by  Cramer 
and  Hansen  [Refs.  2,12].  Each  array  contained  1024  silicon 
photodiodes  on  a  single  chip  with  25  micron  spacing.  The 
accompanying  circuits  provided  a  sampled  and  held  output 
which  was  essentially  analog  except  for  switching  transients. 

At  the  end  of  each  diode  scan  there  was  a  delay  before  the  next 
scan.  During  this  delay  the  diodes  were  reset  and  allowed 
to  measure  the  intensity  of  the  scattering  profile  again. 

The  scanning  of  the  diode  array  repeated  continuously.  The 


•/as  about  34  milliseconds 


actual  sampling  time  of  the  array  v 
with  a  delay  between  scans  of  about  6  milliseconds. 

The  50  centimeter  focul  length  of  the  objective  lens 
combined  with  the  dimensions  of  the  diode  array  provided  a 
half  angle  field  of  view  of  about  3  degrees  for  mediums  of 
refractive  index  near  one.  The  effective  field  of  view  was 
reduced  to  about  2.3  degrees  for  calibrations  when  the 
refractive  index  of  a  Plexiglas  container  and  water  was 
taken  into  account. 

The  laser  beam  collimator  mentioned  previously  at 
first  produced  a  beam  one  centimeter  in  diameter.  A  lens 
in  the  collimator  was  changed  to  reduce  the  beam  diameter  for 
several  reasons.  Extraneous  light  would  be  generated  if 
a  large  beam  impinged  on  the  aperture  of  the  motor  test 
section  window.  Also,  if  the  aft  beam  was  larger  than  the 
motor  exhaust  jet  it  would  be  refracted  in  the  density 
(and  refractive  index)  gradient  between  the  exhaust  and  air. 

The  last  part  of  the  apparatus  was  the  rocket  motor 
itself.  It  was  the  one  used  by  Cramer  and  Hansen  [Refs.  2,12 
and  in  the  present  study  served  only  for  aligning  the  optics. 


IV.  DATA  ACQUISITION'  SYSTEM 


A.  NEW  CONTROLLER 

Hansen  (Ref.  12]  describes  the  major  components  of  the 
Hewlett  Packard  3054A  data  acquisition  system.  A  list  of  the 
manuals  relevant  to  this  study  is  in  Table  1.  The  HP  85 
computer  used  by  Hansen  and  Cramer  was  replaced  with  an 
HP  9836S  as  system  controller.  This  newer  computer  has  far 
more  capability  than  the  HP  85,  including  a  choice  of  more 
powerful  operating  systems.  The  system  used  for  this 
study  was  Basic  Extended  2.1. 

The  data  acquisition  program  written  by  Hansen  needed 
minor  modification  to  acquire  multiple  consecutive  scans  of 
the  photodiodes.  Some  different  I/O  commands  such  as  those 
which  transfer  data  to  the  disk  were  also  incorporated. 

The  revised  version  of  this  program  is  listed  in  Appendix  A. 

A  general  flow  chart  is  presented  in  Figure  3. 

The  9836S  has  two  internal  disk  drives  which  were  used 
to  store  the  data  after  acquisition.  The  data  from  both 
diode  arrays  was  stored  in  the  same  file.  The  eight  (8) 
scans  of  the  motor  cavity  were  first,  followed  by  the  four  (4) 


exhaust  scans. 


3.  MODIFICATIONS 

The  memory  capacity  of  the  .Multiprogrammer  unit  was 
increased  so  that  multiple  consecutive  scans  of  the  diode 
arrays  could  be  recorded  during  a  motor  firing.  This  would 
provide  a  more  statistically  valid  measurement  of  particle 
size.  Fluctuations  of  scattered  light  intensity  for  a 
polydispersion  need  to  be  integrated  over  time  or  averaged 
to  provide  a  more  appropriate  measurement. 

In  order  to  fully  use  the  memory  added  and  make  data 
management  easier  the  data  acquisition  system  needed  to  be 
modified.  The  memory  consisted  of  three  (3)  cards,  each 
with  a  capacity  of  4096  values.  The  fact  that  this  was  a 
multiple  of  1024  (the  exact  number  of  photodiodes)  meant 
that  the  idle  period  between  scans  needed  to  be  excluded 
from  the  data.  If  this  was  not  done,  one  (1)  less  scan  per 
card  would  have  been  acquired  and  locating  the  scans  in 
the  overall  block  of  memory  would  have  been  more  difficult. 

It  was  also  necessary  to  chain  two  of  the  memory  cards 
together  in  a  way  which  would  allow  one  card  to  be  filled 
and  then  the  other.  A  schematic  of  the  data  acquisition 
system  is  presented  in  Figure  4. 

C.  PACING  AND  MEMORY  CONTROL  CIRCUITS 


The  timing  clock  and  blanking  pulse  of  the  photodiode 
circuitry  provided  the  means  for  pacing  data  acquisition. 


Specific  results  of  the  nodi f leathers  were: 

(1)  Memory  space  was  fully  utilized  and  management  of  the 
multiple  scans  made  easier. 

(2)  A/D  conversions  of  the  data  were  made  exactly  when  a 
diode's  output  was  on  line  and  steady.  Thus,  the  analog 
filter  used  in  the  previous  study  was  no  longer 

needed  to  suppress  the  switching  spikes  on  the  data  line. 

The  following  is  a  description  of  the  signals  and  circuits 
used  to  modify  the  data  acquisition  system.  All  voltage 
levels  were  TTL .  A  timing  diagram  in  Figure  5  shows  the 
relations  between  signals.  A  schematic  of  the  circuit  is 
presented  in  Figure  6 . 

The  clock  pulse  was  a  positive  going  spike  at  a  frequency 
of  about  30  KHz.  This  clock  controlled  all  circuits  of 
the  photodiodes.  It  ran  continuously,  even  during  the  blank 
period  between  scans  when  the  diode  output  was  clamped  at 
zero  volts. 

The  blanking  pulse  was  a  signal  which  fell  to  zero  at 
the  beginning  of  each  scan.  It  then  went  positive  at  the 
end  of  the  scan  and  remained  high  until  the  next  scan  began. 

The  clock  pulse  was  used  to  drive  a  pulse  shaper  (mono¬ 
stable  multivibrator) .  This  ensured  that  the  voltage  levels 
through  the  rest  of  the  circuit  would  not  accidentally  fall 
below  the  TTL  threshold.  The  pulse  width  of  the  shaper 
was  adjusted  so  that  the  negative  going  edge  of  each  pulse 
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would  occur  after  the  switching  transient  on  the  data  line 
had  decayed.  This  negative  going  edge  would  eventually 
trigger  the  A/D  converter  to  store  the  output  of  each  diode. 

The  blanking  pulse  was  inverted  and  connected  to  an  AND 
gate  along  with  the  pulse  shaper  output.  The  output  of  this 
AND  gate  is  shown  in  Figure  5  as  the  pulse  shaper  signal  held 
low  between  scans  of  the  diode  array.  This  was  the  basic 
signal  which  paced  data  acquisition. 

This  basic  trigger  signal  was  connected  to  an  AND  gate 
along  with  the  output  of  the  Multiprogrammers  Timer  Pacer 
card.  In  this  way  the  trigger  would  not  reach  the  A/D  until 
the  Timer  Pacer  output  a  pulse.  This  enabling  pulse  from 
the  Timer  Pacer  was  at  least  as  long  as  the  time  for  eight 
scans  of  a  diode  array. 

The  controller  programed  the  Timer  Pacer  to  produce  the 
pulse  when  the  Timer  Pacer  received  a  trigger  from  the  blanking 
pulse.  In  this  way,  data  acquisition  began  at  the  start  of 
a  scan  and  no  data  was  taken  during  the  time  between  scans. 

The  circuit  to  chain  the  memory  cards  together  was 
basically  an  OPv  gate  used  as  a  negative  logic  AND  gate.  The 
end  of  conversion  (EOC)  signal  of  the  A/D  and  the  (FULL)  signal 
of  one  memory  card  were  connected  to  the  gate.  When  both 
signals  went  low  the  second  card  was  then  able  to  store  data 
from  the  A/D.  This  arrangement  is  shown  in  the  schematic  of 
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the  data  acquisition  system  in  Figure  4.  The  automatic  lock- 
out  feature  of  the  memory  cards  when  full,  and  the  relatively 
slow  rate  of  data  throughput  made  it  unnecessary  to  control 
other  handshake  lines  [memory  card  manual,  Table  1]. 

The  circuit  was  designed  to  handle  four  (4)  memory  cards, 
so  no  modification  will  be  necessary  if  one  more  card  is 
added  to  the  system.  This  would  provide  an  additional  four 


scans  of  the  exhaust  beam. 


V.  DATA  REDUCTION 


The  data  reduction  programs  written  by  Hansen  were  not 
used  for  this  study.  The  new  computer  lended  itself  to 
another  approach.  Its  memory  capacity  made  it  unnecessary  to 
chain  programs  together  and  polynomial  curve  fitting  was 
eliminated  in  favor  of  interactive  graphics.  Avoiding 
polynomial  fits  preserved  the  nature  of  the  raw  data  so  that 
one  had  a  better  feel  for  the  parameters.  The  data  reduction 
program  " RDC "  is  listed  in  Appendix  A.  Figure  7  is  a  general 
flow  chart  for  the  program. 

The  following  is  a  description  of  the  program.  The  user 
was  first  prompted  for  values  needed  to  analyze  a  given 
data  set.  For  example,  the  wavelength  of  the  laser  used  and 
the  index  of  refraction  of  the  medium  must  be  known  for 
any  data  set.  Next,  one  had  the  choice  of  reducing  raw 
data  scans  or  reviewing  a  reduced  data  file.  For  raw 
data  one  chose  either  the  exhaust  or  motor  cavity 
beam  data. 

Raw  data  was  plotted  on  the  CRT  and  any  obviously  erroneous 
scans  were  excluded  from  further  reduction.  The  valid  scans 
were  averaged  to  obtain  a  mean  scattering  profile.  The  mean 
intensity  profile  taken  before  particles  were  introduced 
was  then  subtracted  from  that  taken  with  particles  present. 


This  corrected  for  the  characteristics  of  individual  photo¬ 
diodes  and  extraneous  light  which  was  independent  of  the 
particles . 

A  symmetric  moving-average-type  of  digital  filter  was 
then  applied  to  the  profile  to  achieve  some  smoothing. 

This  type  of  digital  filter  was  chosen  for  simplicity  and 
because  it  does  not  have  the  phase  lag  of  analog  filters 
[Ref.  16] .  Preserving  the  phase  of  the  data  was  necessary 
to  retain  angular  resolution.  Another  advantage  of  filtering 
in  the  software  rather  than  hardware  was  that  raw  data 
files  remained  unmolested. 

The  scattering  profile  was  then  analyzed  using  interactive 
graphics.  If  earlier,  one  chosen  to  review  a  reduced  file, 
program  execution  began  here. 

One  had  to  normalize  a  scattering  profile  in  order  to 
compare  it  to  the  theoretical  curves  for  polydispersions. 

The  scattered  intensity  on  the  centerline  of  the  beam  was 

the  correct  value  to  use  for  normalization  but  was  unmeasurable 

due  to  the  beam's  presence. 

The  other  unknown  was,  of  course,  the  particle  size.  These 
two  variables  (centerline  intensity  for  the  measured  profile 
and  D^2  for  the  theoretical  profile  of  normalized  intensity 
vs.  (8^))  were  adjusted  using  interactive  graphics  until 
the  curve  for  polydispersions  coincided  with  the  data.  In 
this  way  the  mean  diameter  of  particles  was  determined. 


The  second  reduction  technique  used  was  the  direct 
application  of  the  method  presented  by  Buchele  [Ref.  13]. 

The  equation  for  the  polydispersion  curve: 

I n ( 0 )  =  EXP  -(.57  aS)2 

was  applied  at  two  points  of  the  scattering  profile.  This 
gave: 

12/11  =  EXP  -D2[(82  -  62)(.57  tt / X )  2 ] 

Solving  this  for  the  diameter  gave: 

D  =  [-Ln(l2/Il)  (X/.  57*)  2/ (62  -  e2)]*5 

The  computer  would  sweep  through  the  data  using  many 
values  of  thetal  (0^  along  with  several  angle  ratios  to 
determine  theta2  ( 0 2 )  .  The  results  were  presented  graphically 
as  particle  size  vs.  thetal  for  each  angle  ratio  (Q Q ^)  . 

In  actual  practice  the  range  of  useable  angles  depends  on 
the  appratus,  and  the  quality  of  the  data.  Therefore,  in  order 
to  interpret  the  results  one  must  have  previously  inspected 
the  data.  The  interactive  graphics  routine  was  well  suited 
to  this  and  provided  a  hard  copy  for  inspection. 

After  reducing  a  set  of  raw  data  the  mean  scattering 
profile  was  stored  on  disk  for  later  review. 


VI.  CALIBRATION  AND  EVALUATION 


A.  IMPROVEMENTS 

The  geometry  of  the  apparatus  used  in  the  investigation 
by  Cramer  and  Hansen  is  compared  with  that  of  this  study 
in  Figure  8.  In  the  previous  study  the  transmitted  beam 
was  allowed  to  enter  the  receiving  optics.  The  beam  was 
focused  off  the  diode  array  a  few  millimeters  from  the  first 
diode.  This  was  necessary  to  avoid  damaging  the  diodes 
but  introduced  some  uncertainty  in  angle  measurements. 

The  intense  image  of  the  beam  along  with  scattered  light  from 
the  receiving  optics  produced  a  high  level  of  extraneous 
light.  In  the  present  study,  stops  were  used  to  intercept 
the  beam  before  reaching  the  receiving  optics.  These  stops 
provided  several  advantages.  A  high  intensity  beam  could 
be  used  while  producing  little  extraneous  light.  Optics 
alignment  was  also  improved.  This  reduced  error  in  angle 
measurement.  Alignment  was  accomplished  using  a  neutral 
density  filter  to  reduce  beam  intensity  and  protect  the 
diodes.  A  schematic  of  the  apparatus  is  in  Figure  2.  The 
laser,  collimator,  beam  splitter  and  prism  were  positioned  so 
that  the  beams  passed  through  the  appropriate  measurement 
areas.  The  narrow  pass  filter  and  imaging  lens  were  then 
positioned  so  that  the  beam  entered  on  the  centerline.  The 


photodiode  array  was  then  moved  using  a  three-axis  micrometer 
so  that  the  focused  beam  fell  on  the  first  diode.  The  beam 
stops  were  then  put  in  place  and  the  neutral  density  filter 
removed.  In  this  way  measurements  commenced  exactly  from 
the  optical  axis  of  the  beam. 

Procedures  were  also  refined  to  account  for  the  bending 
of  light  rays  as  they  passed  through  the  walls  of  the  particle 
container.  As  noted  earlier,  the  index  of  refraction  of  the 
container  and  the  medium  containing  particles  affects 
scattering  measurements.  A  Plexiglas  box  held  the  particle 
samples  and  a  magnetic  stirrer  kept  the  samples  suspended 
in  water.  The  index  of  refraction  of  the  Plexiglas  and 
water  combination  was  measured  using  a  simple  technique. 

A  microscope  was  used  to  measure  the  ratio  of  actual  depth 
to  apparent  depth  for  Plexiglas  and  water.  The  index  of 
refraction  was  determined  to  be  1.39.  This  value  was 
applied  to  the  data  to  convert  the  measured  scattering  profile 
to  that  actually  produced  in  the  medium  containing  the 
particles . 

B.  RESULTS 

Calibration  results  are  summarized  in  Table  2.  Initial 
tests  were  done  with  two  samples  wihch  were  basically  mono¬ 
dispersions  of  large  particles.  Figure  9  shows  the  measured 
profile  of  scattered  light  for  glass  spheres  ranging  from 


37  to  44  microns  in  diameter.  This  profile  was  obtained  by 
placing  the  focussed  beam  just  far  enough  from  the  first 
diode  to  avoid  saturation  with  no  particles  present.  The 
diodes  located  at  angles  less  than  about  .01  radians 
saturated.  The  first  bright  ring  for  particles  of  about  40 
micron  diameter  was  visible  near  .02  radians.  Figure  10 
shows  a  profile  for  the  same  particles,  illustrating  use  of 
the  beam  stops  to  avoid  diode  saturation  and  improve  angle 
measurements.  In  this  case,  the  first  diode  was  located 
exactly  on  the  centerline  of  the  beam  as  discussed  above. 

Results  for  a  sample  of  53  to  63  micron  glass  spheres  are 
shown  in  Figure  11.  The  center  lobe  was  nearly  completely 
missed  but  the  first  two  bright  rings  were  seen  near  .014 
and  .022  radians.  The  first  two  dark  rings  near  .01  and 
.019  radians  were  also  seen.  The  method  described  earlier 
of  setting  the  beam  spread  parameter  (a9)  equal  to  the 
zeroes  and  maximums  of  J^(a9)  was  used  to  calculate  a  size 
of  about  58  microns.  Also  shown  in  Figure  11  is  the 
theoretical  profile  for  a  polydispersion  with  =  54 

microns . 

Various  polydispersions  of  either  glass  spheres  or 


aluminum  oxide  powder  were  then  tested.  These  polydispersions 
consisted  of  fairly  large  particles.  Results  are  shown  in 
Figures  12  through  17.  These  tests  showed  that  the  apparatus 


had  two  distinct  modes  of  operation.  If  the  particle  concen¬ 
tration  was  very  high,  or  if  large  particles  dominated  the 
polydispersion,  many  of  the  diodes  at  the  smaller  angles 
would  saturate.  This  left  only  the  data  at  larger  angles 
useable.  When  many  diodes  saturated,  the  theoretical  curve 
given  by  Dobbins,  et  al.,  was  used  to  determine  size.  This 
was  done  because  this  curve  was  valid  for  the  larger  angles 
and  lower  relative  intensities.  The  curve  from  Buchele 
[Ref.  13]  was  not  valid  for  values  of  the  beam  spread 
parameter  greater  than  three  (3)  . 

For  low  particle  concentrations  and/or  small  particles  the 
data  proved  more  accurate  at  the  smaller  angles.  If  no 
diodes  were  seen  to  saturate  then  one  knew  the  measurement 
was  in  the  higher  intensity  part  of  the  center  lobe.  Here 
the  curve  given  by  Buchele  was  quite  satisfactory  for  sizing. 

The  smallest  particles  tested  were  five,  ten,  and  twenty 
micron  polystyrene  spheres.  The  bright  rings  for  these 
particles  occurred  at  angles  too  large  for  the  apparatus  to 
measure.  For  these  samples  the  diodes  did  not  saturate. 

Both  the  Gaussian  curve  fit  and  the  two  angle  method  were 
used  to  obtain  032*  These  results  were  especially  consistent. 
It  should  be  noted  that  the  two-angle  method  uses  the  equation 
for  the  Gaussian.  If  the  measured  profile  matched  the 
Gaussian  exactly,  then  would  be  the  same  for  any  (6^) 

and  angle  ratio  employed.  Some  variations  in 


calculated  -ue  to  the  imperfect  fit  are  obvious  ir. 

Figures  19,  21  and  23. 

A  scanning  electron  microscope  was  used  to  photograph 
the  types  of  particles  tested.  These  photographs  are  shown 
in  Figures  24  through  27.  Equation  (1)  was  used  along  with 
these  photographs  to  calculate  some  of  the  values  of 
in  Table  II.  Calculations  of  D32  for  the  polystyrene  were 
arrived  at  using  the  manufacturers  data  on  size  distributions 
The  photographs  generally  confirm  the  validity  of  the 
technique . 


VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  results  of  the  calibration  tests  showed  that  the 
apparatus  is  capable  of  accurately  measuring  mean  particle 
size  for  a  broad  range  of  mean  diameters.  It  was  found  that 
the  technique  was  most  accurate  if  the  theoretical  profile 
fit  or  the  two-angle  method  were  applied  at  the  smallest 
possible  scattering  angles. 

The  rocket  exhaust  is  likely  to  attenuate  the  beam 
somewhat,  reducing  the  problems  related  to  diode  saturation 
at  small  scattering  angles.  Thus,  measurements  should  be 
possible  using  the  high  intensity  part  of  the  center  lobe. 
This  should  make  data  reduction  less  ambiguous.  Actual 
testing  should  begin  with  measurements  at  the  exhaust 
plane  of  the  motor.  These  should  be  compared  with  collected 
exhaust  samples  to  validate  the  use  of  the  apparatus  in 
an  actual  motor  environment.  Measurements  in  the  motor 
cavity  would  then  be  interpreted  based  on  the  correlation 
between  exhaust  samples  and  exhaust  measurements. 

It  is  also  recommended  that  the  index  of  refraction  of 
the  combustion  gases  be  investigated.  A  literature  search 
for  an  estimate  of  the  index  of  refraction  would  probably 
be  satisfactory. 


TABLE  I 


ELECTRONICS  MANUALS 

1.  HP  Memory  Cards  Model  6970B  Operating  Manual 

2.  HP  Timer  Pacer  Card  Model  69737A  Operating  Manual 

3.  HP  Analog  to  Digital  Converter  Card  Model  69736A 
Operating  Manual 

4.  HP  Users  Guide,  "Using  the  9826  and  9836  Computers  with 
the  6942A  Multiprogrammer" 

5.  Basic  Language  Reference  Guide  with  Extensions  2.0  for 
Series  200  Computers 
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TABLE  II.  CALIBRATION  RESULTS 


Figure  2.  Schematic  of  I.ight  Scattering  Apparatus. 


Figure 


Figure  5.  Timing  Diagram  for  Data  Acquisition . 


Schematic  of  Pacing  Circuit  for  Data  Acquisition 


Figure  7.  (Continued)  Flow  Chart  for  Program  "RDC 
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Figure  7.  Flow  Chart  for  Program  "RDC" . 


Geometries  for  Light  Scattering  Measurements 
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Micron  Glass  Spheres  Using  No  Beam  Stop 


CURVE  FIT  RESULTS 
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Micron  Glass  Spheres  Using  Beam  Stop 


CURVE  FIT  RESULTS 
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Micron  Glass  Spheres  Showinq  Diffract 


CURVE  FIT  RESULTS 
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37  Micron  Glass  Spheres,  Hiqh  Concentra 


CURVE  FIT  RESULTS 
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Micron  Aluminum  Oxide,  Low  Concent 
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cron  Aluminum  Oxide,  Hiqh  Concentration 


CURVE  FIT  RESULTS 
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Figure  17.  25  Micron  Aluminum  Oxide,  Low  Concentration. 


CURVE  FIT  RESULTS 
INTENSITY  vs.THETR 


Figure  18.  5  Micron  Polystyrene,  Curve  Fit. 


RNGLE  METHOD 


Figure  19.  5  Micron  Polystyrene,  Two  Angle  Method 
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Figure  20.  10  Micron  Polystyrene,  Curve  Fit. 


ANGLE  METHOD 


Figure  21.  10  Micron  Polystyrene,  Two  Angle  Method 


CURVE  FIT  RESULTS 
INTENSITY  vs.THETR 


Figure  22.  20  Micron  Polystyrene,  Curve  Fit. 
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Figure  23.  20  Micron  Polystyrene,  Two  Angle  Method 
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APPENDIX  A 
PROGRAM  LISTINGS 

ihhhhh  ACQDTA  :  FOR  9S38S  ****************** 

liiiHiim  ACQUIRES  MULTIPLE  SCAMS  OF  BOTH  DIODE  ARRAYS  «* »«******«»** 
iihhhhi  AND  STORES  THEN  ON  DISK  ***************** 

IMHHHH  BERT  HANSEN,  KEN  GRAHAM  .KELLY  HARRIS  1984  *t**m********** 

OPTION  BASE  1 

COM  91*1121.92*112]  ,Tli[28I  ,Dl*tl2I  .D2SI 12]  ,Address<4),E(4096) 

COM  Ad 024)  BUFFER, B14096)  BUFFER ,C(4Q9fc)  BUFFER  ,D(409S)  BUFER 
OUTPUT  709;'AR"  ! ANALOG  RESET 

ASSIGN  SMultilO  TO  72310  ! CLEARS  THE  WAKE-UP  SERVICE  REQUEST 

ENTER  gMult il 0 ;Qql  ,9q2,9q3)9a4,9c5,Qqti  !CF  THE  HULTIPRQGRAMMER 

HAT  Address=  (0) 

ENABLE  INTR  7 
CLEAR  722 
PRINT  USING  "3/* 

PRINT  "ENTER  THE  FILENAMES  OF  THE  DATA  FILES  TO  EE  CREATED  (e.o.  RAW1,RAW2 

PRINT  "  12  CHARACTERS  MAXIMUM,  EACH’ 

PRINT  USING  7" 

PRINT  "  AN  EMPTY  DISK  MUST  BE  ^  THE  LEFT  DISK  DRIVE  • 

PRINT  •" 

PRINT  "  THE  DATA  FILES  WILL  NEARLY  FILL  A  DISK" 

PRINT  USING  7" 

MASS  STORAGE  IS  INTERNAL ,4,0"  !  CHANGE  THIS  LATER  IF  NECESSARY 

INPUT  "  INPUT  FILE  NAMES  NOW  -  (FILENAME! ,FILENAME2>  ",D1*,D2* 

M5=1024 

Zz$=" : INTERNAL ,4 , 1  *  !  STRING  INDICATES  MASS  STORAGE 

CREATE  BOAT  DlHZz*,6144,16  !  12  SCANS  OF  1024  9  2*9  BYTES  PER  RECORD 
CREATE  BDAT  D2*4Zz*, 6144,16  !  1024*12/2=NUMBER  OF  REC0RDS=6144> 

PRINT  USING  "?" 

PRINT  *  DATA  HILL  BE  STORED  ON  DISKETTE  HITH  FOLLOWING  FILE  NAMES:" 
PRINT  USING  "///" 

PRINT  "  NO  PARTICLES - FILENAME  =  ";D1* 

PRINT  "  PARTICLES - FILENAME  =  ";D2* 

PRINT  USING  "//" 

PRINT  "  IS  THIS  A  CALIBRATION  ?  ENTER  '  Y  '  IF  YES  " 

PRINT  "  ANYTHING  ELSE  IF  NO  " 

INPUT  R* 

PRINT  USING  "//" 

PRINT  "  BE  SURE  LASER  IS  ON" 

PRINT  “  PRESS  [CONTINUE]  WHEN  READY" 

PAUSE 


•  •  •.  ■;  -.-y  ■-*  ’ 


VV 


410  !*=■  y  then  ro 

423  PRINT  USING  '3' 

438  INPUT  'ENTER  THE  THRESHOLD  PRESSURE  TO  TRIGGER  THE  DEVICES  <psi)\S8 

440  INPUT  'ENTER  TINE  DELAY  FROM  THRESHOLD  PRESSURE  (sec) ',T8 

458  !  TIMES  ARE  INTERPRETED  BY  THE  COMPUTER  IN  SECONDS  DOWN  TO  ,001 

460  V0=S3/151.5  !  CALIBRATION  OF  PRESSURE  TRANSDUCER 

470  D3i=Dlt  !  NO-PARTICLES  STRING  NAME 

488  GOSUB  Multiprog 

498  GOSUB  Storedata 

500  IF  R$='Y*  THEN  1000 

510  !************************************************************************ 

520  Contcheck:  ICont inoit y  checking 

530  OUTPUT  709;'AC20'  ! CONNECT  3456A  DVM  TO  IGNITION  FIRING  LINE 

540  LOCAL  LOCXOUT  7  IKEEPS  OPR,  FROM  SWITCHING  FRONT  PANEL  TO  LOCAL 

550  R7=225  !  IN  THOUSANDS  OF  OHMS 

560  OUTPUT  722; 'HSM0O2SW2SO! LI S0F4R 1 T3QX I '  ! RESISTANCE  MEASUREMENT 

570  !  H  RESETS  DVM;SM802  SETS  SERVICE  REQUEST  MASK  WHERE  002  IS  OCTAL  REP,  OF 
580  !  THE  SERIAL  POLL  MASX  BYTE;SU2  TELLS  WHICH  TERMINAL  SWITCH  IS  USEDjSOl 

590  !  SYSTEM  OUTPUT  MODE  ON  -  WAITS  FOR  CONTROLLER  TO  HANDSHAKE ;L1  LOAD 

600  !  INTERNAL  MEMORY  ON;SO-FUNCTION  SHIFT  CFF-.F4-TW0  WIRE  CONNECTION  TO  DUN 
610  !  R 1 -AUTORANGING ; T3-SINGLE  TRIGGER ;Q-LOAD  INTERNAL  MEMORY  OFF ;X1 -EXECUTE 

620  !  STORED  PROGRAM. 

630  GOSUB  Reading 
640  RS=V 

650  IF  R8<R7  THEN  GOTO  Ok 

660  ! *******************»**CONTINUITY  CHECK********************************* 

670  PRINT  USING  *8' 

680  PRINT  'CONTINUITY  CHECK  BAD!!!!' 

690  PRINT  'RECHECK  BEFORE  PROCEEDING.  WHEN  CHECKED,  PRESS  [CONTINUE I* 

700  BEEP  3008 ,.3 

710  WAIT  .1 

720  BEEP  100,1.0 

730  PAUSE 

740  GOTO  Contcheck 

750  OkiPRINT  USING  '0' 

760  DISP  '  CONTINUITY  CHECK  IS  O.K.' 

770  PRINT  '  BE  SURE  NITROGEN  IS  ON' 

780  OUTPUT  709;'AC21'  iCONNECT  DVM  FOR  IGNITOR  VOLTAGE  MEASUREMENT 

790  PRINT  '  DVM  CONNECTED  TO  FIRE  SWITCH' 

800  PRINT  USING  V 

810  PRINT  'BE  SURE  VISICORDER  IS  SET  UP  TO  RUN  ON  PROPER  SCALE  WITH  LAMP  ON.' 
820  PRINT  USING  V 

830  OUTPUT  722 ; 'HSMO 02SW2Z0SO1L1FL0 , 0 1 STIS0F1 R 1 T3Q '  !  VOLTAGE  ON  IGNITER 

840  DISP  '  STANDING  BY  FOR  IGNITION' 

850  PRINT  '  STANDING  BY  FOR  IGNITION' 

860  BEEP  2000,. 1 

870  OUTPUT  722; *X1'  ! VOLTMETER  TRIGGER 

880  GOSUB  Reading  !READS  VOLTMETER 

890  R9=ABS(V) 

901  IF  R9U0  THEN  GOTO  870  !  12  VOLTS  ON  IGNITOR 


!  VOLTAGE  ON  IGNITER 


(VOLTMETER  TRIGGER 
!READS  VOLTMETER 


!  12  VOLTS  ON  IGNITOR 


*  *  *  ■  •'*  «"■  **•  ,*• 

_•/ V a- o V  oV  .  * 
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910 

920 

933 

940 

958 

961 

970 

980 

998 


C0=TlfE’A!E 
OUTPUT  709; 'AC22* 
OUTPUT  722;*Xr 
C0SU6  leading 
R9=ABS(V) 

IF  R9< VO  THEN  GOTO  930 
WAIT  T8 
01=T1HEDATE 
IF  MOT  THEN  1060 


!  USc.S 


COMPUTER  CISC*  TO  SET  ELAPSES  TINES 
!  CONNECT  PRESSURE  XDUCER  TO  DVM 
!  TRIGGER  VOLTMETER 


!  THRESHOLD  PRESSURE 
!  TINE  DELAY 


! IF  ACTUAL  RUN  THEN  SKIP  SOME  LINES 
!  FOR  CALIBRATIONS  INTRODUCE  PARTICLES  AND  THEN  CONTINUE 

PRINT  CHRIC12) 

PRINT  USING  ■////////• 

PRINT  *  INTRODUCE  PARTICLES  THEN  PRESS  KEY  *  9  TO  TAKE  DATA’ 

ON  KEY  9  LABEL  ’  TAKE  DATA*  GOTO  1060 
1050  Standby:  GOTO  Standby 
1060  GOSUB  Multiprog 


1010 

1028 

1030 

1040 


!  TAKE  PARTICLE  DATA 


1070 

1080 

1090 

1100 

1110 

1120 

1130 


! PARTICLE  DATA  FILE  NAME 


ELAPSED  TIME: 
ELAPSED  TIME: 


FIRE  TO  MULTIPROGRAMMER  CALL  =  *;01-00 
FIRE  TO  MEMORY  INTERRUPT  =  ,:02-00 


DATA  STORED  ON  DISK  WITH  FILENAMES  (’;D1$;*)  AND  ( “ ■, D25 ; * ) * 


LOCAL  7 
D3$=D2$ 

GOSUB  Storedata 
PRINT  * 

PRINT  * 

PRINT  ■ 

GOTO  End 
1148  Multiprog;  ! 

1150  OUTPUT  723; ‘00,2,3, 12, 13T* 

OUTPUT  723;,CC,5,6,9,10T*  ! CLEARS  THE  ARM, BUSY  AND  ECP  OF  MEMORY  CARDS 
OUTPUT  723; *CC, 1,4,1 IT’  ! CLEARS  SAME  FOR  A  TO  D'S  AND  TIMER  PACER 
OUTPUT  723;*SF,2,3)1,.001 ,12,T'  !  THE  (1)  IS  2'S  COMPLIMENT  BINARY 

OUTPUT  723;*SF,5,3,1,.001,12,T*  !  THE  .001  IS  THE  LEAST  SIGNIFICANT  BIT 

OUTPUT  723;,SF,9,3,1,.001,12,T1  !  THE  12  IS  FOR  12  BIT  WORD  SIZE 

OUTPUT  723;‘SF,12,3,1 , .00t,12,T*  !  SINCE  THE  A  TO  D  IS  12  BIT 
OUTPUT  723;-WF,3,1023T,WF,6,4095T’ 

OUTPUT  723 ; -UF , 1 0 , 4095T , WF , 1 3 , 4095T ’ 

!  SETS  REFERENCE  WORD  FOR  WHEN  TO  STOP  TAKING  DATA  AND  GENERATE  INTERRUPT 
OUTPUT  723^,2. 1, IT, WF,5.1,1T,WF, 9.1, IT, WF, 12.1, IT1  ‘SETS  FIFO  MODE 
OUTPUT  723i-WF,3.1,OT,WF,3.2,OT,WF,3.3,OT,WF,13.1,OT,WF,13.2,OT,WF,13.3,OT 


1160 

1170 

1180 

1190 

1280 

1210 

1220 

1230 

1240 

1250 

1260 


1270  OUTPUT  723;'WF,6.1,OT,WF,6.2,0T,WF,6.3,OT,WF,10.1,0T,WF,10.2,0T,WF,10.3,0T 


1280  !  SETS  COUNTER  AND  POINTERS  OF  2ND  MEMORY  CARD  IN  EACH  PAIR  TO  0 
1290  OUTPUT  723;'AC,3T,AC,6T,AC,10T,AC,I3T'!ARMS  CARDS  WHICH  GENERATE  INTERRUPT 
1300  OUTPUT  723 ; "WF , 4 .2, OT ,WF , 4 , 1ST"  !  TIMER  PACER  GIVES  1  PULSE  OF  1  SEC  WHEN 
1310  !  TRIGGERED  BY  THE  BLANKING  PULSE  (PLENTY  OF  TIME  FOR  8  SCANS) 
1328  WAIT  2 

1330  K=SPOLL(723)  IWAIT  FOR  MEM  INTERRUPT 
1340  IF  K 064  THEN  GOTO  1330 

1350  OUTPUT  723;,WF,4.2,1T*  !  MAY  BE  UNNECESSARY  TO  ALTER  TIMER  PACER  SINCE 
1360  !  MEMORY  CARDS  HAVE  AUTOMATIC  LOCKOUT  BUT  FOR  NOW  WE  WILL  DO  IT 
1370  02=TIMEDATE  ' 

1380  ON  ERROR  GOTO  Err.trap  !  NEEDED  TO  READ  ARMED  CARD  INTERRUPT  LIST 
1390  SEND  7;UNL  NLA  TALK  23  SEC  12  !  SPECIFICALLY  ASKS  FOR  INTERRUPT  LIST 


*.  -v 


1433  Vir _rejd :  ESTER  7, Address!* 
1413  flencaras:  PR IHT  1 
1420  PRINT  • 

1430  NAT  Address*  (0) 

1440  OFF  ERROR 

1450  OUTPUT  723;,DC,3,6,13,10T* 
1460  OUTPUT  723) 'HI ,2,1 024T ■ 

1470  ENTER  72305; A<*> 

1480  OUTPUT  723; ,MI,5J4096T* 

1490  ENTER  72305;B(») 

1500  PRINT  * 

1510  OUTPUT  723 ;*MI, 12,40961' 

1520  ENTER  72305;C(») 

1530  OUTPUT  723;,HI,9,4096T* 

1540  ENTER  72305;D(*> 

1550  ENABLE  INTR  7; 8 

1560  PRINT  *  I 

1570  MAT  A=  (-1)*A  !  THE 

1580  HAT  B=  (-1  )«B  !  PERFl 

1590  MAT  C=  (-1)»C 

1600  HAT  D=  (-1)*D  !  DIODI 

1610  RETURN 

1620  Storedata : ! 

1630  ASSIGN  PDiskf lie  TO  D3*iZz* 
1640  ASSIGN  BBuffl  TO  BUFFER  C<») 
1650  ASSIGN  RBwff2  TO  BUFFER  D<*> 
1660  ASSIGN  @Buff3  TO  BUFFER  B(») 
1670  CONTROL  PBuffl, 3)1,32763,1 
1680  CONTROL  BB«ff2,3;l ,32767,1 
1690  CONTROL  @B«ff3, 3)1,32767,1 
1700  TRANSFER  RBoffl  TO  SDiskfile 
1710  WAIT  FOR  EOT  BDiskfile 
1720  TRANSFER  8B»ff2  TO  BDiskfile 
1730  WAIT  FOR  EOT  BDiskfile 
1740  TRANSFER  8Buff3  TO  8Diskfile 
1750  WAIT  FOR  EOT  BD1Skfile 
1760  ASSIGN  BDiskfile  TO  * 

1770  ASSIGN  BBuffl  TO  * 

1780  ASSIGN  BBuff2  TO  * 

1790  ASSIGN  BBiff3  TO  * 

1800  RETURN 

1810  Reading:! 

1820  STATUS  7, ljAO 
1830  ENTER  722;V 
1840  ENABLE  INTR  7;8 


i  (READ  WHICH  CARDS  INTERRUPTED 

MEMORY  CARDS  WHICH  GENERATED  INTERRUPTS  ARE  * 
SLOTS*  =  *;Address(») 


IDISSARH  HEM  CARDS 
!  SET  UP  CARD  TO  BE  READ 
!  GETS  DATA  FROM  1024  MEMORY  FOARD 


EXHAUST 


!  GETS  DATA  FROH  4096  MEMORY  BOARD 
DATA  ENTERED* 


MOTOR  DATA  ENTERED* 

!  THE  1024  CARD  IS  INCLUDED  BUT  NOT  SAVED.  IT  DIDN'T 
!  PERFORM  WELL.  COULD  BE  REPLACED  BY  A  4096  CARD. 

!  DIODE  VOLTAGES  ARE  NEGATIVE  SO  SIGNS  ARE  CHANGED 


!  MOTOR  CAVITY  4  SCANS 
!  MOTOR  CAVITY  4  SCANS 
!  EXHAUST  4  SCANS 
‘SETS  BUFFER  POINTERS  TO  FULL 
! INTERFACE  REGISTERS  SECTION  OF 
{LANGUAGE  MANUAL 
!  ORDER  OF  DATA  ON  THE  DISK  IS 
!  MOTOR  CAVITY --8  SCANS 
!  EXHAUST  -  4  SCANS 
!WAIT  BECAUSE  OVERLAPPING 
! TRANSFERS  ARE  NOT  WANTED 

! CLOSE  I  /  O  PATHS 


! CHECKING  STATUS  BEFORE  READING 
'VOLTMETER  IS  A  FORMALITY  TO 
! CLEAR  THE  SERVICE  REQUEST 


1850  RETURN 

1860  Err.trap:  IF  ERRN=159  AND  ERRUVar  read)  THEN  Hemrards 

1870  PRINT  ERRH*  SEVEN  IF  THE  ERROR  WAS  NOT  THE  ONE  PLANNED 

1880  GOTO  Mencards  !FOR  PROGRAM  EXECUTION  CONTINUES 

1890  End;  END 
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1 0  1  mHHHitHHHiilllHHHKtHHiHHHHHHHtHHIHHHHIH 

20  liHIHimiHHHtl  ROC  XXXXXXXXXXXXXXXXXXXX 

30  Ixxxxxxxxxxxxxxxxxxx  PLOTS  RAW  DATA  xxxxxxxxxxxxxxxxxxxx 
40  !x«xxxxxxxxxxxxxxxxx  AVERAGES  xxxxxxxxxxxxxxxxxxxx 

50  !xxxtxxxxxxx«xxxxxxx  FILTERS  xxxxxxxxxxxxxxxxxxxx 

60  ixxxxxxxxxxxxxxxxxxx  DETERMINES  MEAN  DIAMETER  xxxxxxxhxixxxxxxxxx 

70  ixxxxxxxxxxxxxxxxxxx  BY  INTERACTIVE  GRAPHICS  xxxxxxxxxxxxxxxxxxxx 

90  Ixxxxxxxxxxxxxxxxxxx  AND  THE  TWO-ANGLE  METHOD  xxxxxxxxxxxxxxxxxxxx 

90  ixxxxxxxxxxxxxxxxxxx  Robert  Kell v  Harris  xxxxxxxxxxxxxxxxxxxx 

100  Ixxxxxxxxxxxxxxxxxxx  1984  xxxxxxxxxxxxxxxxxxxx 

111  | XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

120  OPTION  BASE  I 

130  COM  /Hrdaaus/  Av2(1024) 

140  COM  /Gauss/  T1 ( 1 024) ,G( 1024) ,L 

151  COH  /Mai/  H7,H5,Xt,Yl,Xfl,YH 

160  COH  /Readata/  BlPlH,Q3$t201,04M21]>ZiM20],Yl(8192)  BUFFER, Y2( 81 92)  BUFFE 

R 

170  COM  /Two/  Avl (1024> ,M,M1 ,F 

180  DIM  5cans(S),X(1024) 

191  COM  /Plots/  P1I120] ,P2$[20]lP3$t201,P4iI201 

200  INTEGER  Graf (1:12480)  BUFFER 

210  Choose:  PRINT  CHRIU2) 


220 

Old=0 

230 

PRINT  ’  TO  LOOK  AT  ANY  DATA  FILE  THE  PROGRAM  NEEDS  SOME  STARTING  INFORHATI 

ON* 

240 

PRINT  " 

250 

PRINT  * 

TO  ACCOUNT  FOR  THE  CHANGE  IN  WAVELENGTH  IN  THE  MEDIUM’ 

260 

PRINT  ” 

270 

PRINT  1 

ENTER  THE  INDEX  OF  REFRACTION  OF  THE  MEDIUM’ 

280 

PRINT  " 

290 

PRINT  * 

WATER  =1.33’ 

300 

PRINT  ' 

AIR  =  1.0’ 

310 

PRINT  * 

ESTIMATE  OF  EXHAUST  =  1.1’ 

320 

BEEP 

330 

INPUT  " 

THIS  VALUE  ADJUSTS  THE  COMPARISON  CURVE  TO  THE  MEDIUM’ ,M1 

340 

PRINT  CHRM12) 

351 

PRINT  ’ 

TO  ACCOUNT  FOR  REFRACTION  OF  LIGHT  AT  BOUNDARIES  BETWEEN’ 

360 

PRINT  *• 

370 

PRINT  1 

THE  MEDIUM  AND  AIR  YOU  MUST  ENTER  THE  INDEX  OF  REFRACTION  * 

380 

PRINT  •* 

C?  THE  COMBINATION  OF  THE  MEDIUM  AND  ITS  BOUNDARY 


390 

PRINT  *  C?  THE  COMBINATION  OF 

THE  MEDIUM  AND  ITS  BOUNDARY1 

400 

PRINT  “ 

410 

PRINT  *  THIS  VALUE  IS  APPLIED  DIRECTLY  TO  THE  DATA* 

420 

PRINT  ** 

430 

PRINT  *  AIR=1.0 

PLEXIGLASS  6  WATER  =  1.39* 

440 

PRINT  '• 

450 

PRINT  ‘ESTIMATE  OF  MOTOR  CAVITY  COMBUSTION  PRODUCTS  6  WINDOW  =  1.22* 

460 

PRINT  " 

470 

PRINT  ‘  ESTIMATE  OF  ROCKET  EXHAUST  =  1.1  OR  1.0‘ 

460 

INPUT  M 

490 

PRINT  CHRK12) 

500 

INPUT  ‘  ENTER  LASER  WAVELENGTH  <HeNe=. 6328. Ar=. 488)* ,L 

510 

INPUT  ‘ENTER  THE  DESIRED  PLOTTING  INTERVAL  OF  DIODE  ARRAY  (2,4,6)‘,H 

520 

F=460 

Inn  FOCUL  LENGTH  OF  OBJECTIVE  LENS 

530 

Diod-.025 

Inn  DIODE  SPACING  MAY  BE  .03 

540 

FOR  1=0  TO  1023 

550 

X(IM)=I+1 

! CREATE  AN  ARRAY  OF  DIODE  NUMBERS 

560 

Tl(I+l)=(I«Diod)/(F*M) 

! CREATE  AN  ARRAY  OF  THETA 

570 

NEXT  I 

560 

PRINT  CHR$(12) 

590 

CALL  Displayl(Old) 

I  JUST  PRINTS  INITIAL  REMARKS  ON  CRT 

600 

IF  Old  THEN  CALL  Review! Avl (*)) 

!  THESE  TWO  LINES  ALLOW  FOR 

611 

IF  Old  THEN  GOTO  Gauss 

! REVIEWING  REDUCED  DATA  FILES 

620 

BEP 

630 

INPUT  *  INPUT  TWO  FILENAMES  NOW  (NO-PART  ,PART)‘,Q3$,Q4$ 

640 

INPUT  ‘  ENTER  M'  FOR  MOTOR 

'2'  FOR  EXHAUST ‘,P 

650 

ZjO=,:INTERNAL,4)1‘  ! 

USES  LEFT  DISK  DRIVE  FOR  RAW  DATA  FILES 

660 

CALL  Readata 

670 

CALL  Shift 

6S0 

BEEP 

690 

Screen:  ! 

700 

PRINT  CHR$(12) 

(CLEARS  SCREEN 

710 

PI $= ‘VOLTAGE* 

720 

P2$=*DI0DE  NUMBER' 

!  STRINGS  FOR  PLOTS 

730 

P3t= ‘VOLTAGE  vs.  DIODE* 

740 

P 40= ‘RAW  DATA' 

750 

M5=1024 

760 

Xt=64 

!  SET  UP  VALUES  FOR  PLOTS 

770 

Tt=.04 

780 

Xn=4 

790 

Yn=5 

800 

CALL  Plot (t2( ») ) 

!  DRAWS  AXES  ,ETC... 

810 

GSTORE  Craf(*> 

!  STORES  GRAPHICS  DISLAY  JUST  MADE 

820 

CALL  Dataplot(B,Yl(*) ,H) 

'PLOT  NO-PARTICLES  DATA 

830 

PRINT  ‘ 

NO-PARTICLES  PLOTTED’ 

840 

PRINT  " 

S8S  PRINT  *  WHICH  SCANS  ARE  GOOD  1. . .ie. . .1 ,2,4, 5,7, INCLUDE  LAST  CO 

m‘  .'LAST  COrtHA  IS  REQUIRED  OR  YOU  HAVE  TO  HIT  CONT. .  TWICE 

890  INPUT  Scans!*) 

900  PRINT  CHR4U2) 

91ft  PRINT  USING  •/////////' 

920  PRINT  '  AVERAGING  THE  SELECTED  SCANS' 

930  GRAPHICS  OFF 

940  CALL  AveraaetJ, Scans!*), Yl!*),Avl(*)) 

950  PRINT  USING 

961  GLOAD  Graf!*)  ILOADS  GRAPHICS  ARRAY  RATHER  THAN  HASTE  TIME  RE-DRAUING 
970  GRAPHICS  ON 

980  CALL  Result (Aul (*),X!*),H)  !  PLOTS  1024  ELEMENT  ARRAYS 

990  PRINT  USING  7/' 

1000  PRINT  1  Average  Intensity  No-Particles' 

1010  ON  KEY  0  LABEL  ' AVER AGE 1'  GOTO  Screen 
1020  ON  KEY  1  LABEL  'PLOT-PARTICLES*  GOTO  1060 

1030  PRINT  1  PRESS  KEY  *  0  TO  RE-AVERAGE  OR  *  1  TO  CONTINUE* 

1040  Standby:  !  HANUAL  CALLS  THIS  INTERRUPT  DRIVEN  PROGRAMMING 

1050  GOTO  Standby  !  LOOPS,  WAITING  FOR  USER  TO  DECIDE 

1060  PRINT  USING  '8* 

1070  OFF  KEY  1  ! HELPS  AVOID  CONFUSION  BY  CLEARING  THAT  BOX 

1080  CALL  Dataplot(B,Y2(*),H)  'PLOT  PARTICLES  DATA 

1090  PRINT  USING  V 

1100  PRINT  '  PARTICLE  DATA  PLOTTED' 

1110  PRINT  *  FOR  A  HARD  COPY  OF  THIS  RAH  DATA  PRESS  KEY  *  6" 

1120  PRINT  '  ' 

1130  PRINT  1  OR  TO  CONTINUE  PRESS  KEY  »  1' 

1140  ON  KEY  6  LABEL  'HARD  6  RAH*  GOTO  Raw 
1150  ON  XEY  1  LABEL  'CONTINUE  '  GOTO  Select 
1160  GOTO  Standby 

1170  Raw:  CALL  Plot ( Y2( *) f 1 )  !  THE  ONE  (1)  IS  AN  OPTIONAL  PARAMETER 

1180  CALL  Dataplot(B,Y2!*),H)  !  WHICH  IS  USED  TO  GET  HARD  COPIES 

1190  Select:! 

1200  OFF  KEY  6 

1210  OFF  KEY  1 

1220  PRINT  CHR0U2) 

1230  PRINT  USING  '////' 

1240  PRINT  '  HOW  MANY  SCANS  SEEN  TO  BE  GOOD  ?' 

1250  INPUT  J 

1260  PRINT  '  WHICH  SCANS  ARE  GOOD  ?...ie... 1,2, 4,5, 7, INCLUDE 

LAST  COHHA' 

1270  INPUT  Scans!*) 

1280  PRINT  CHRIU2) 

1290  PRINT  USING  '/////////' 

1300  GRAPHICS  OFF 

1310  PRINT  '  AVERAGING  THE  SELECTED  SCANS' 

1320  CALL  Average! J, Scans!*) ,Y2!») ,Av2<») ) 

1330  PRINT  USING  '8* 

1340  P4*='AVERAGED  SCANS' 


1353 

CALL  Plott Av2(») > 

1340 

GRAPHICS  ON 

1370 

GSTORE  Graft*) 

1380 

CALL  Res»lttAv2!*),X(*),H) 

1390 

PRINT  USING  •//' 

1400 

PRINT  ' 

Averane  Scattered  Intensity* 

1410 

ON  KEY  2  LABEL  ’AVERACE2* 

GOTO  1080 

1420 

ON  XEY  3  LABEL  'SUBTRACT' 

GOTO  1450 

1430 

PRINT  ' 

HIT  XEY  0  2  TO  RE-AVERAGE  OR  KEY  *  3  TO 

CONTINUE' 

1440 

GOTO  Standby 

1450 

PRINT  USING  '?' 

1440 

OFF  XEY  3 

1470 

HAT  Aul=  Av2-Avl 

{SUBTRACTS  NO-PARTICLES  FROM  PARTICLES 

1480 

GLOAD  Graft*) 

1490 

CALL  Result(Av1(*),X(*),H) 

1500 

PRINT  USING  '///' 

1510 

PRINT  ' 

Plit  of  the  Difference  Between  Particles  and  N 

•-Particles' 

1520 

PRINT  ' 

XEY  ♦  4  FOR  HARD  COPY* 

1530 

PRINT  ' 

XEY  0  1  TO  CONTINUE  * 

1540 

ON  XEY  4  LABEL  '  HARD  AVERAGE'  GOTO  1570 

1550 

ON  XEY  1  LABEL  '  FILTER 

'  GOTO  1400 

1540 

GOTO  Standby 

1570 

CALL  Pl»t(A»l (*), 1 ) 

1580 

CALL  Result (Avl (*) ,X( *) ,H) 

1590 

PLOTTER  IS  3, 'INTERNAL* 

1400 

PRINT  USING  '//' 

1410 

PRINT  * 

ENTER  THE  NUMBER  OF  TINES  YOU  UISH  TO  APPLY* 

1420 

PRINT  ' 

TO  DIGITAL  FILTER  FOR  SMOOTHING* 

1430 

PRINT  ' 

EXAMPLE********  10* 

1440 

PRINT  ' 

TAXES  ABOUT  1.5  MINUTES* 

1450 

INPUT  ' 

ZERO  IS  ALSO  ACCEPTABLE *,Fil 

1440 

IF  Fil=0  THEN  Ga»ss 

1470 

P4t= 'FILTERED  DATA* 

1480 

CALL  FiltertA«l(*),Fil) 

1490 

PRINT  CHR$<12) 

1700 

CALL  Piet (Avl (*)) 

1710 

IF  Old  THEN  HAT  X=  T1 

1720 

CALL  Reseltt Avl (»),X(*),H) 

1730 

PRINT  ' 

Piet  of  the  Difference  Between  Particles  and  N 

•-Particles* 

1740 

TER' 

PRINT  ' 

AFTER  APPLICATION  OF  A  DIGITAL  FIL 

1750 

PRINT  * 

for  a  HARD  COPY  * 

1740 

PRINT  ' 

PREPARE  the  PLOTTER  and* 

1770 

PRINT  ' 

PRESS  KEY  0  8* 

1780 

PRINT  ' 

OR* 

1790 

PRINT  * 

PRESS  KEY  »  9  TO  NORMALIZE* 

1800 

ON  XEY  8  LABEL  'HARD  FILTERED*  GOTO  1830 

181*  OH  KEY  9  LABEL  *  NORMALIZE  *  GOTO  Gavss 


!  The  1  is  for  hard  copy 


1820  GOTO  Stancbv 
1830  CALL  Plut(Avl(») .!) 

1840  CALL  Res»lt(Awl<*M(*),H) 

1850  OFF  KEY  8 

1880  Gauss:  OFF  KEY  2  ! KEY  1  OPTION  IS  STILL  VALID 

1870  OFF  KEY  8 

1880  ON  KEY  1  LABEL  1  FILTER  '  GOTO  1600  !  ALLOWS  ONE  TO  SMOOTH  OLD  DATA 
1890  !  WHEN  REVIEWING  REDUCED  FILES 

1900  PRINT  CHR$( 12) 

1918  !A«2  Amy  is  Normalized  BUT  Avl(«)  is  Still  Saved  For  Re-work  if  Needed 
1920  HAT  Av2=  Avl  !  TO  BEGIN  WITH  THE  ARRAY  IS  ASSUMED  TO  BE  NORMALIZED 
1930  !  AND  THE  USER  CHANGES  THIS  WITH  INTERACTIVE  GRAPHICS 

1940  !  IN  THE  SUBROUTINE  "Conoare' 

1950  Pl$= 'NORMALIZED  INTENSITY* 

1960  P2t=*THETA  (rad)'  ISTRINGS  FOR  PLOTS 

1970  P3*='INTENSITY  us. THETA" 

1980  P4$= 'CURVE  FIT  RESULTS' 

1990  H5=.05  (SET  UP  VALUFS 

2000  Xt=,004  IFOR  PLOTS 

2010  Yt=.l 
2020  X«=2 
2030  Y«=2 

2040  CALL  Plot(Av2<»>> 

2050  CALL  Result(Av2(*),Tl<*>,H)  (PLOT  OF  NORMALIZED  INTENSITY  us,  THETA 

2060  GSTORE  GraflO  !  SAVES  THE  SCREEN  IMAGE 

2070  CALL  Compare! Avl (OjHjMjMl, Graf <*)) 

2080  GRAPHICS  OFF 

2091  ON  KEY  4  LABEL  'OTHER  ARRAY'  GOTO  Choose 
2100  ON  KEY  5  LABEL  '  TWO-ANGLE"  GOTO  Bochele 

2110  ON  KEY  9  LABEL  *  NORMALIZE'  GOTO  Gauss 

2120  ON  KEY  6  LABEL  'STORE  DATA'  GOTO  Storedata 

2130  ON  KEY  7  LABEL  '  QUIT  *  GOTO  Quit 

2140  CALL  Menuel  !  PRINTS  OPTIONS  ON  THE  SCREEN 

2150  GOTO  Standby 
2160  Buchele:  ! 

2170  CALL  Twoanqle 

2180  CALL  Henuel 

2190  GOTO  Standby 

2200  Storedata:! 

2211  CALL  Store 

2220  CALL  Henuel 

2230  GOTO  Standby 

2240  Quit:  ! 

2250  END 

2260  SUB  Aweraoe! J,Scans(*) ,Y(*> ,Av(*)>  'AVERAGES  SELECTED  SCANS 
2270  HAT  Av=  (0)  UNITIALIZES  THE  ARRAY  LOCAL  TO  THIS  ROUTINE 

2280  FOR  1=1  TO  J  STEP  1 

2290  K=(Scans(I)-l 1*1024+1  !  THIS  COUNTER  IS  THE  BEGINNING 

2300  L=Scans(I)»1024  !  AND  THIS  ONE  THE  END  OF  BLOCKS 

2310  FOR  I1=K  TO  L  !  OF  1024  IN  THE  OVERALL  DATA 


!  THIS  COUNTER  IS  .ALWAYS  BETWEEN 
!  I  AND  1024 


2320 

I2=IH(+1 

!  THIS  COUNTER  IS  ALWAYS  BETWEEN 

2339 

A»(I2)=4#(I2)+Y(I1)  !  1  AND  1024 

2340 

NEXT  11 

2350 

NEXT  I 

2369 

HAT  A*=  A#/(J) 

!*****THE  AVERAGE*****! 

237* 

SUBEND 

2380 

SUB  Dataplot(B,Y(*) ,H) 

>9  IS  4  OR  8 (THE  NUMBER  OF  SCANS) 

2399 

COH  /Max/  H7,H5,Xt,Yl,Xn,Yii 

2410 

LDIR  9 

2411 

LORG  4 

2420 

FOR  J=1  TO  B 

!EACH  SCAN 

2430 

MOVE  509+3*50, M7-. 95 

!  SEE  NEXT  LINE 

2441 

LABEL  J;  •  HELPS  KEEP  TRACK  OF  EACH  SCAN  AS  IT  APPEARS 

2459 

LINE  TYPE  1 

2469 

K=(J-1)*1024+1 

(BEGINNING  OF  EACH  SCAN  AND 

2479 

L=J*1924 

!THE  END  WITHIN  THE  TOTAL  BLOCK 

2489 

MOVE  1,Y(K) 

!  HOVE  TO  THE  FIRST  POINT 

2490 

FOR  I=K  TO  L  STEP  H 

2500 

X=I-K+1 

!  THIS  GIVES  1  TO  1924  FOR  ABSCISSA 

2510 

PLOT  X,Y(I) 

2520 

NEXT  I 

2530 

NEXT  J 

2549 

PENUP 

2559 

SUBEND 

2560 

SUB  Result(Iav(*),D(*),H) 

•THE  AVERAGE  INTENSITY  IS  PLOTTED 

2570 

clip  ON 

2580 

MOVE  D(l),Iav(l) 

2590 

FOR  1=1  TO  1024  STEP 

H 

2600 

DRAW  D(l),Iav(I) 

2610 

NEXT  I 

2620 

PENUP 

2630 

SUBEND 

2640 

SUB  Plot(Y(»), OPTIONAL  Hard) 

•IF  THE  OPTIDNAL(Hard) IS  RECEIVED 

2650 

COM  /Max/  M7,M5,Xf ,Yt,Xrt,Y« 

!THE  FIGURE  GOES  TO  THE  PLOTTER 

2669 

COH  /Plats/  Plt[20),P2M20],P3*[2fl],P4$l20] 

2670 

M7=MAX(Y(*))*1,1 

<  A  SCALING  VARIABLE 

2689 

M6=.1*M7 

!  ANOTHER  SCALING  VARIABLE 

2690 

PRINT  " 

2700 

GINIT 

IGINIT  IS  JUST  GOOD  PRACTICE  SO  YOU  KNOW 

2710 

INHERE  YOU  ARE  BEGINNING 

2720 

SELECT  NPAR 

•THIS  DETECTS  IF  THE  HARD  COPY  IS  DESIRED 

2730 

CASE  1 

! IT  COULD  BE  DONE  WITH  IF  THEN  LOGIC 

2740 

PLOTTER  IS  3, ’INTERNAL* 

! BUT  IS  PRESENTED  FOR  FAMILIARIZATION 

2750 

CASE  2 

•SINCE  IT  IS  MORE  POWERFUL  IN  COMPLEX 

2760 

PLOTTER  IS  795,‘HPGL' 

•SITUATIONS 

2770  END  SELECT 

2780  GRAPHICS  OH 

!  TO  BE  ABLE  TO  SEE  THE  PLOT 

2790  CSI2E  5,. 6 

!  TO  BE  ABLE  TO  SEE  THE  PLOT 


!  SCREEN  UNITS  FOR  KARGINS 
!  DRAWS  A  BOX 
!  EXTENT  CF  X  AND  Y 


2820  LORG  5 
2830  MOVE  75,95 
2840  LABEL  USING  'X';P4* 

2850  MOVE  75,90 
2860  LABEL  P3$ 

2870  MOVE  75,20 
2880  LABEL  P2* 

2891  LDIR  90 
2900  HOVE  18, B0 
2910  LABEL  PI* 

2920  VIEWPORT  30,125,30,85 
2930  FRAME 

2940  WINDOW  0,N5,-H6,M7 
2950  AXES  Xt,M6,0,0,X*,l,5 
2940  LDIR  0 
2970  CSI2E  3,, 5 
2980  LORG  8 
2991  CLIP  OFF 

3010  FOR  I=-N4  TO  M7  STEP  M7/10  INUMBER  THE  Y  AXIS 

3010  MOVE  0,1 

3020  LABEL  USING  *0,HD.DDDD';I 

3030  NEXT  I 

3040  CSIZE  3,. 6 

3050  LDIR  90 

3040  LORG  8 

3070  FOR  1=0  TO  M5  STEP  Xt  INUMBER  THE  X  AXIS 

3080  MOVE  I, -HA 

3090  LABEL  USING  ,*,X';I 

3100  NEXT  I 

3110  PENUP 

3120  SUBEND 

3130  SUB  Conpare(Avl(*),H,M, Ml, INTEGER  Graft*)) 

3140  COM  /Gauss/  T1(*),G(*),L 

3150  COM  /Hrdgaus/  Av2<») 

3140  COM  /Max/  H7,H5,Xt,Yt,X«,Y« 

3170  DIM  T(l:1024),AtU:17),Dobbins(l:l7),Tld(l:17> 

3180  DATA  3. ,3.5,4, ,4.5,5. ,5. 5, 4., 4. 5, 7, ,7,5,8, ,8, 5, 9,, 9. 5,11.  ,10.5,1 1 . 

3190  DATA  .0554,. 035,. 0204,. 014,. 0104,, 0081,. 00605,. 005, .00374,. 0032,. 00248, 

.0022, .00185, .0014, .00135, .0012, .001 
3200  READ  AtU) 

3210  READ  Dobbins! *) 

3220  OFF  XEY  (GETS  RID  OF  ALL  LABELS  ON  KEYS 

3230  ON  KEY  3  LABEL  'MENUE'  GOTO  Sobexit 

3240  PLOTTER  IS  3, 'INTERNAL'  !  IN  CASE  A  HARD  COPY  WAS  JUST  MADE 

3250  H*rd=0  !  (0)  SO  ONE  DOESN'T  EXIT  TOO  SOON 

3240  CenterlineM  I  THE  INITIAL  NORMALIZING  VALUE 

3270  D=20.  !  INITIAL  PATICLE  MEAN  DIAMETER  IN  MICRONS 

3280  Change:  CSIZE  4,. 4 
3290  PRINT  USING  '//////’ 


3200 

POINT  * 

STAND  BY  FOR  CURVE* 

3310 

HOVE  0,1 

3320 

CLIP  ON 

3330 

Con=PI*D*Hl/L 

!A  CONSTANT 

3340 

FOR  1=1  TO  1024  STEP  H 

3350 

T(I)=Tl(I)*Con 

I  Theta-bar  for  stne  Dbarlnote  HI) 

3360 

G(I)=EXP<-(( .57*T(I))*2)) 

IGaussian  curve  for  this  Dbar 

3370 

IF  T(I)>3.0  THEN  3400 

■GAUSSIAN  NO  GOOD  FOR  Theta-bar >3 

3380 

DRAW  T1(I),G(I) 

IGaussian  curve  vs.  theta 

3390 

NEXT  I 

3400 

HAT  Tld=  At/tCon) 

! THETA  FOR  POINTS  ON  DOBBINS 'CUR 

VE 

3410 

FOR  1=1  TO  17 

3420 

DRAU  Tld < I) .Dobbins! I ) 

3430 

NEXT  I 

3440 

PENUP 

3450 

IF  Hard  THEN  SUBEXIT  IEXITS  ROUTINE  IF  A  HARD  COPY  WAS  JUST  HADE 

3460 

PRINT  USING  -/////////■ 

3470 

PRINT  * 

USE  THE  KNOB  TO  VARY  THE  PARTICLE 

SIZE* 

3480 

PRINT  * 

OR  HIT  KEY  »  6  FOR  HARD  COPY* 

3490 

PRINT  *  * 

3500 

PRINT  ' 

KEY  *  9  ALLOWS  YOU  TO* 

3510 

PRINT  ■ 

RE-NOR  HALIZE* 

3521 

PRINT  *  ' 

3530 

PRINT  * 

HIT  KEY  »  3  TO  GET  OUT* 

3540 

PRINT  USING  '//////■ 

3550 

ON  KEY  6  LABEL  ’HARD  COPY*  GOTO  Hardgauss 

3560 

ON  KEY  9  LABEL  'NGRHALIZE*  GOTO  Divide 

3570 

ON  KNOB  .05  GOTO  Pulse 

H.85)  IS  TIHE  INTERVAL  IN  WHICH 

3580  Wait:  GOTO  Uait 

!  PULSES  FROH  THE  KNOB  ARE 

3590  Pulse:  POINT  USING  *8* 

!  COUNTED  AND  THIS  NUHBER  IS 

3600 

Strngl$=*Hean  Diaaeter  * 

!  USED  BY  THE  INTERRACTIVE 

3610 

Strng2$="  Hicrins* 

!  GRAPHICS  TO  VARY  THE  PARTICLE 

3620 

C»unt=KNOBX 

!  SIZE  AND  PLOT  THE  ASSOCIATED 

3630 

D=DROUND ( D+C«  un  t/1 5 , 2) 

!  GAUSSIAN  APPROXIMATIONS  OF 

3640 

GLOAD  Graft*) 

!  SCATTERING  PROFILES 

3650 

LDIR  0 

3660 

LORG  8 

3670 

HOVE  H5,.8*H7 

3680 

LABEL  StrnglUVALt(D)4Strng2* 

3690 

GOTO  Change 

3700  Hardoauss: ! 

3710 

PRINT  * 

PREPARE  the  PLOTTER* 

3720 

PRINT  ■ 

PRESS  CONTINUE  fir* 

3730 

PRINT  * 

HARD  COPY* 

3740 

PAUSE 

3750 

Hard=l 

!  SO  THAT  SUBEXIT  OCCURS  AFTER  HARD  COPY 

3760 

CALL  Plot <Av2( *) , 1 ) 

!  1»<HARD  COPY 

3770 

CALL  Result<Av2(*),Tl(*),H) 

3783 

PRINT  USING  ‘8* 

3790 

LORG  8 

3801 

MOVE  M5,.8»H7 

3810 

LDIR  9 

3820 

LABEL  StrngmVAL»(D)iStrng2* 

3830 

HOVE  0,1 

3840 

GOTO  Change 

3850  Divide ; l 

3860 

PEN  -1 

(ERASES  LAST  PLOT 

3870 

HOVE  0,0 

3880 

FOR  1=1  TO  1024  STEP  H 

3898 

DRAW  T1(I),Av2(I) 

3900 

NEXT  I 

3910 

PRINT  CHR9U2) 

3920 

PRINT  USING  '//////■ 

3930 

PRINT  1 

USE  THE  KNOB  TO  RE-NORHALIZE* 

3940 

PRINT  ' 

KEY  *  4  WILL  RE-PLOT  USING1 

3950 

PRINT  * 

UPDATED  SCALE1 

3960 

ON  KEY  4  LABEL  1  RE-PLOT1 

GOTO  Replot 

3970 

ON  KNOB  .05  GOTO  Vary 

3980 

GOTO  Wait 

3990  Vary: PRINT  CHR$<12) 

4000 

PEN  -1 

(ERASES  LAST  PLOT 

4010 

HOVE  0,1 

4021 

FOR  1=1  TO  1024  STEP  H 

4031 

ORAN  Tl(I),Av2(I) 

4041 

NEXT  I 

4051 

PEN  1 

4060 

C2=KN0BX 

4070 

Centerline=Centerline+C2/120 

(KNOB  USED  TO  VARY  NORMALIZING  VALUE 

4080 

NAT  Av2=  Avl/lCenferline) 

(NEU  NORMALIZED  DATA  ARRAY 

4090 

GOTO  3870 

(PLOTS  THE  NEU  NORMALIZED  ARRAY 

4100  Replot:  ! 

4111 

PRINT  CHR$(12) 

4120 

CALL  Piet (Av2(* ) ) 

4130 

CALL  Result(Av2(*),Tl(*),H) 

4140 

GSTORE  Graft*) 

4150 

GOTO  Change 

4160  Subexit:  ! 

4170 

SUBEND 

4180 

SUB  Henvel 

4190 

PRINT  USING  'V 

4200 

PRINT  1 

YOU  CAN  RE-NORHALIZE1 

4218 

PRINT  'imilliHtmiimHIiiHIKtKKtimiHIHOHHHHHIlt1 

4220 

PRINT  1 

YOU  CAN  RE-AVERAGE  (New  Data  Only)1 

4230 

PRINT  'HIHIHIIIIHtHHIHimiHHHHitlllllHIHHHtHIHmil 

4240 

PRINT  * 

PRESS  KEY  0  4  TO  LOOK  AT  OTHER  DATA1 

4250 

PRINT  * 

OLD  OR  NEW  1 

4261 

PRINT  1 

MOTOR  /  EXHAUST1 

4271 

PRINT  'IHHiHHiHIIHIHHIIIIHtiimtIHIHHIHIHIHiiHXli 

4533 

PRINT  • 

PRESS  KEY  #  5  ’ 

4293 

PRINT  1 

FOR  THE  '  TWO-ANGLE  '  METHOD' 

4390 

PRINT  1 

OF  PARTICLE  SIZING*  * 

4310 

PRINT  'imo******************************************************' 

-  *  — 

4323 

PRINT  ‘ 

TO  STORE  THE  REDUCED  DATA  PRESS  KEY  *  6 

9 

4330 

PRINT  'IHHHHHOiiHHiiHKHiHKOHHIKMHOiHOtiOtHH' 

4340 

PRINT  * 

TO  END  THIS  SESSION  PRESS  KEY  t  7' 

4350 

SUBEND 

4380 

SUB  Displayl (Old) 

:v  - 

4370 

PRINT  1 

TO  REDUCE  NEW  DATA  PRESS  KEY  *  1' 

-  -  ■ :  : 

4380 

PRINT  " 

1 

4390 

PRINT  ’ 

TO  REVIEW  PREVIOUSLY  REDUCED  DATA  PRESS  KEY  *  2* 

4400 

PRINT  *  • 

4410 

ON  KEY  1  LABEL  * 

NEW  DATA'  GOTO  New 

4420 

ON  KEY  2  LABEL  ‘ 

OLD  DATA'  GOTO  Review 

4430  Wait:  GOTO  Wait 

4440  New:  PRINT  CHRS412 

) 

1 

4450 

PRINT  • 

PUT  THE  DISK  IN  THE  LEFT  DRIVE  IF  IT  IS  NOT  ALREADY' 

4480 

PRINT  " 

4470 

PRINT  * 

ENTER  THE  NAMES  OF  THE  FIRST  AND  SECOND  FILES.' 

4480 

PRINT  " 

4490 

PRINT  1 

EACH  FILE  HAS  DATA  FROM  BOTH  DIODE  ARRAYS.' 

4500 

PRINT  " 

9 

4510 

PRINT  * 

FIRST  IS  NO-PARTICLES  6  NEXT  IS  PARTICLES  - D1$,D2*' 

4520 

PRINT  ** 

v.;;' 

4530 

Old=0 

4541 

SUBEXIT 

S-S-S- 

4550  Review;  01d=l 

! THIS  VARIABLE  IS  PASSED  TO  THE  MAIN  PROGRAM  TO  INDICATE 

* 

4580 

ITHAT  THE  DATA  TO  BE  READ  IS  ONE  (1)  SCAN  AND  THAT 

w 

<  *v 

4571 

ITHE  AVERAGING  ROUTINES  WILL  NOT  BE  USED 

»  *  . 

.  ‘'v' 

4580 

SUBEND 

4590 

SUB  Readata 

•  • 

.*  **  •' 

4800 

COH  /Readata/  B,P,H,Q3$l28],Q4$l201,Zz$[201,Yl(»)  BUFFER, Y2(»)  BUFFER 

4810 

Xyz=l 

!  FILE  POINTER  VARIABLE 

i 

4620 

B=8 

!  NUMBER  OF  SCANS  IN  MOTOR  DATA 

p 

4630 

R 1=65536 

INUMBER  OF  BYTES  OF  HOTOR  DATA 

4640 

IF  P=2  THEN  Xyz=4097  ! RECORD  0  WHERE  EXHAUST  DATA  BEGINS 

4650 

IF  P=2  THEN  R 1=32768  INUMBER  OF  BYTES  OF  EXHAUST  DATA 

‘S'. 

4663  IF  P=2  THEM  8=4  ! NUMBER  OF  SCAMS  IN  EXHAUST  DATA 

4870  PRINT  *  " 

4880  PRINT  *  READING  DATA  FROM  FILE  ON  DISK* 


469# 

ASSIGN  RDiskl  TO  Q3$6Zzt 

4730 

ASSIGN  8Disk2  TO  Q4$4Zz» 

IOPEN  I/O  PATHS 

4710 

ASSIGN  BBgffl  TO  BUFFER  Y!<*) 

4723 

ASSIGN  ?B«ff2  TO  BUFFER  Y2<*> 

4730 

CONTROL  0Diskl,5;Xyz 

4740 

TRANSFER  RDzskl  TO  RBuffl ;COUNT  R1 

IREADS  NO-PARTICLE  DATA 

4750 

WAIT  FOR  EOT  SDiskl 

I  NOT  AN  OVERLAPPING  TRANSFER 

4763 

CONTROL  8Disk2,5;Xyz 

IXyz  SETS  DISK  FILE  POINTER  TO 

4770 

(EITHER  MOTOR  OR  EXHAUST  DATA 

t 


(READS  PARTICLE  DATA 


4780  TRANSFER  RDisk2  TO  ?Euff2;CCUNT  R1 


4790  WAIT  FOR  EOT  ?Disk2 
4890  ASSIGN  BDiskl  TO  » 

4810  ASSIGN  8Disk2  TO  * 

4820  ASSIGN  SBuffl  TO  * 

4830  ASSIGN  8Buff2  TO  * 

4840  SUBEND 

4850  SUB  Plot2(E,D(»),X> 

4860  DIN  Titlel$[25],ThetjlI203 

4870  GINIT 

4880  IF  X=1  THEN  PLOTTER  IS  705,’HPGL’ 

4890  DEG 

4900  GRAPHICS  ON 

4910  VIEWPORT  35,125,30,85 

4920  Max=10*INT((MAX(D(*))+10)/10) 

4930  WINDOW  DU,2),D(E,2),0,Nai 

4940  G=(D(E,2)-D(1,2))/(E-1)»4 

4950  F=(INT(E/4)-l)*4 

4960  IF  F=0  THEN  F=4 

4970  CLIP  D(4,2)-G,D<F,2)+2*G,0,Hax 

4980  GRID  G,5,D(4,2)-G,I 

4990  CLIP  OFF 

5000  LORG  8 

5010  LDIR  90 

5020  CSIZE  4, .5 

5030  FOR  1=4  TO  E  STEP  4 

5040  MOVE  D<I,2),0 

5050  LABEL  USING  ’.DDDD'jDlI^) 

5060  NEXT  I 

5070  LDIR  0 

5080  LORG  8 

5090  FOR  1=10  TO  Max  STEP  10 

5100  MOVE  D(4,2)-G,I 

5110  LABEL  USING  •*,X'jI 

5120  NEXT  I 

5130  CSIZE  6,. 6 

5140  Title9= ’TWO-ANGLE  METHOD* 

5150  Titlel$='For  Various  Anqle  Ratios 

5160  SizeO=*SIZE  (nicrons)’ 

5170  Thta*=*THETA  (rad)’ 

5180  Sub»=T 

5190  LDIR  90 

5200  LORG  5 

5210  B=D(4,2)-G-(D<E,2)-Dd,2))/10 

5220  MOVE  B, Max/2 

5230  LABEL  Sue* 

5240  LDIR  0 

5250  A=<D(E,2)+Dd,2»/2 

5260  MOVE  A, -Max/4 

5270  LABEL  ThtaO 


(JUST  GOOD  PRACTICE  TO  CLOSE 
(I/O  PATHS 

(DEGREES  for  LABa  DIRECTION 

(AN  X  GRID  LINE  EVERY  4th  POINT 

(MAXES  GRID  UNIFORM 

(PUTS  NUMBERS  ON  X  AXIS 

(NUMBER  Y  AXIS 

(STRINGS  FOR  PLOTS 


5280 

LORG  3 

5290 

MOVE  A, -fU «/4 

5300 

LABEL  Sub* 

5310 

LORG  5 

5320 

MOVE  A, Haul. 15 

5330 

LABEL  Title* 

5340 

CSIZE  4.5 

5350 

MOVE  A, Maul. 05 

5360 

LABEL  Titlel* 

5370 

PENUP 

5381  SUBEND 

5390  SUB  Distribution^*), Tratio,E) 

5400 

LINE  TYPE  1 

5410 

CLIP  ON 

5420 

Tratio*=VAL$(Tr atio ) 

!  ANGLE  RATIO 

5430 

MOVE  D(1 ,2) ,D(1 , 

1) 

5440 

FOR  1=1  TO  £ 

(THIS  SUBROUTINE  PLOTS  THE 

5450 

DRAW  D(I, 2) 

>J>U,1) 

(PARTICLE  SIZES  DERIVED  USING 

5460 

NEXT  I 

(VARIOUS  ANGLE  RATIOS  APPLIED 

5470 

LORG  4 

(TO  THE  DATA  OVER  A  RANGE  OF 

5480 

LINE  TYPE  1 

(ANGLES 

5490 

CSIZE  4,. 3 

5510 

MOVE  DU-1, 2), DU-1, l)+2 

5510 

LABEL  Tratio* 

5520 

PENUP 

5530  SUBEND 

5540  SUB  Twoanqle 

5550 

5560 

!«hihhhhh« 

SUBPROGRAM  UWOANGLE’ 

M»f  S***MX*H**«X* 

5570 

j  **»*«H»«*f  **H 

PARTICLE  SIZING  BY 

5580 

MEASURING  THE  RATIO 

5590 

OF  INTENSITY  AT 

5600 

liHHHHHHH 

TWO  ANGLES 

HiiHHHKiHHH* 

5611 

5620  OPTION  BASE  1 

5630  COM  /Two/  A«1(1024),M,N1,F 

5640  COM  /Gauss/  T1U024),G(1024),L  ‘THE  GAUSS  IS  NOT  USED  HERE 

5650  !  BUT  THE  COM  BLOCK  HAS  THETA 

5660  DIM  0(200,2)  ! AND  WAVELENGTH 

5670  PRINT  CHRIU2) 

5680  PRINT  USING  •///////* 

5690  PRINT  1  THIS  SUBPROGRAM  USES  THE  TWO-ANGLE  METHOD  DESCRIBED  BY  BUCHELE 
5700  PRINT  *  ' 

5710  PRINT  1  TO  CALCULATE  PATICLE  SIZE  FOR  VARIOUS  ANGLE  RATIOS  AND  ANGLES. 
5720  PRINT  *  * 

5730  PRINT  *  IT  IS  HOPED  THAT  THE  CURVES  WHICH  RESULT  WILL  SHOW  ’ 

5740  PRINT  1  ' 

5750  PRINT  *  WHICH  SIZE  IS  THE  MOST  PROBABLE.1 

5760  PRINT  1  1 

5770  PRINT  'AFTER  NOTING  FROM  THE  RAW  DATA  WHICH  ANGLES  CONTAIN  THE  CENTER  LOBE 


TO  CALCULATE  PATICLE  SIZE  FOR  VARIOUS  ANGLE  RATIOS  AND  ANGLES.' 


IT  IS  HOPED  THAT  THE  CURVES  WHICH  RESULT  WILL  SHOW  1 


WHICH  SIZE  IS  THE  MOST  PROBABLE.1 
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5780  PRINT 
5790  PRINT 
5800  PRINT 
5810  PRINT 
5820  PRINT 
5830  PRINT 
5840  INPUT 
5850  01 
5860 


IT  *  ENTER**»**********THE  SMALLEST  USEABLE  ANGLE,' 

IT  "  THE  SMALLEST  ANGLE  RATIO, AND' 

IT  '  THE  STEP  SIZE  BETWEEN  ANGLE  RAIDS' 

IT  *  YOU  WISH  TO  EXPLORE* 

IT  '  EXAMPLE*************  .012,1.2,. 4' 

IT  '  ENTER  thetal ,theta-rJtia,theta-step*,Q,A ,B 

OFF  KEY 

X=0  (GRAPH  WILL  APPEAR  ON  SCREEN 


5870  Begin:  (RUNNING  CONTINUES  HERE  WHEN  A  HARO  COPY  IS  DESIRED 

5880  (Van  De  Hulst  and  Gunprech  6  Sleepevich  explain  that  the  chanae  in 

5890  (wavelength  of  the  bean  is  accounted  for  by  dividing  by  the  index 

5980  *  of  refraction  of  the  nediun.  *********L=L/H1  ********** 

5910  C=  ( L/N 1  / .  57/P  I )  *  2  ! - see  page  15  of  nasa  tech  paper  2156 

5920  FOR  N-l  TO  1024  (to  see  this  is  a  convenient  constant 


5930  IF  Q<T1(N)  THEN  5950 
5940  NEXT  N 

5950  FOR  Tratio=A  TO  3  STEP  B 


5960  FC 
5970 
5980 
5990 
6000 
6018 
6020 
6030 
6040 
6050 
6060 
6070 
6080 
6090 
6100 
6111 
6120 
6138 
6140 
6150 
6160 
6170 
6180 

6190  Xconp: 
6200 


FOR  J=N  TO  1 01 O/Tratio  STEP  10 


(FINDS  POSITION  OF  MINIMUM  ANLGE 

‘VARIOUS  ANGLE  RATIOS 
(SETS  THE  RANGE  OF  POSSIBLE 
!  ANGLES  (THETA  1) 


Th  ta  1  =T1  <  J> 

Thta2=Tl(J*Tratio)  ! 

Deltatheta=Thta2*2-Thtal  *2 
11=0 
12=0 

FOR  I=J-5  TO  J+5  !T 

Il=Il+Avl(I)  (I 

NEXT  I  (O 

11=11/11 

FOR  I=INT(J*Tratio-5)  TO  INT(J*Tratio+5) 


!  THETA  2 


(THESE  2  LOOPS  AVERAGE  A  FEW 
(INTENSITIES  IN  THE  HOPE 
(OF  A  STEADIER  CALCULATION 


I2=I2*Avl(I) 
NEXT  I 

12=12/11 

E=(J-N)/10+1 


!  DONE  HERE  FOR  THETA2  FOR 
!  THE  GIVEN  ANGLE  RATIO 


E=(J-N)/18+1  (THIS  IS  A  COUNTER  FOR  THE  ARRAY  CONTAINING 

(PARTICLE  SIZE  AND  THETA1  FOR  THE  GIVEN 
(ANGLE  RATIO 

IF  II <12  OR  11=0  OR  I2<=0  THEN  D(E,I)=0  (ALLOWANCE  FOR  IF  THE  DATA 
IF  11(12  OR  11=0  OR  I2<=0  THEN  GOTO  Xconp  (IS  NOT  WELL  BEHAVED 
D(E,l)=SOR(-C/Delta the ta*LOG! 12/11 ))  (TWO  ANGLE  METHOD 

(VALUE  OF  DIAMETER  based  on  INTENSITY  RATIO 
(For  a  given  ANGLE  RATIO  and  ANGLE  THETA1 


conp:  D(E,2)=Thtal 

Spar=PI*D(E,l)*Ml/L  (THIS  IS 
Tbar=$par*Thta2 
IF  Tbar>3  THEN  J=1010/Tratio 
IF  E=1  AND  Tbar  >3  THEN  6330 
(the  above  line  ends  all  calculation  if 


(THIS  IS  PI*D/LAMBDA  ~  THE  SIZE  PARAMETER 
(THETA  BAR  FOR  THE  LARGER  ANGLE 
atio  (THIS  ENDS  THE  DO  LOOP  FOR  THIS 

330  (ANGLE  RATIO  SINCE  THE  GAUSSIAN 

ation  if  (IS  NOT  VALID  WHEN  Tbar  )  3 


6258  (the  first  elenent  (E=l)failed  the  test 


FIRST  TIME  THROUGH—  D  HAS  THE  MOST  ELEMENTS  IT  WILL  HAVE 


6230  !  S3  SET  UP  THE  GRAPH  USING  D's  PRESENT  PARAMETERS 

6290  PRINT  CHRtll2) 

6310  IF  Tratio=A  THEN  CALL  Plot2(E,D<») ,X) 

6319  CALL  Distribution(D(*),Tratia,E) 

6320  NEXT  Tratio 


6330  ON  XEY  2  LABEL  ‘HARD  COPY’  GOTO  Hard 

6340  ON  XEY  3  LABEL  ’MENUE  1  GOTO  Sibeiit 

6350  PRINT  USING  ’/////* 

6360  PRINT  ’  YOU  CAN  GET  A  HARD  COPY  BY  PRESSING  XEY  1  2’ 

6370  PRINT  1  OR  EXIT  THIS  ROUTINE  BY  PRESSING  XEY  *  3" 

6389  Standby ‘.GOTO  Standby 

6390  Hard : X=1  !A  VARIABLE  TO  CONTROL  PLOTTER 

6400  GOTO  Begin 

6410  Subexit :  G1NIT 
6420  GRAPHICS  OFF 

6430  SUBEND 
6440  SUB  Shift 

6450  COM  /Readata/  B,P,H, 93*120], 040120], Z:*120],Y1(»)  BUFFER, Y2(»)  BUFFER 
6460  DIN  Ed: 9192) 

6470  PRINT  CHRSC12) 


6480  PRINT  USING  ’///////* 

6490  PRINT  ’  RAW  DATA  IS  BEING  SHIFTED  TO  CORRECT  FOR  SHALL  GAPS  BETWEEN  SCANS’ 
6500  PRINT  *“ 

6510  PRINT  ’  BE  PATIENT,  IT'S  A  LONG  SET  OF  DO-LOOPS’ 

6520  !  THERE  ARE  SOME  SMALL  GAPS  BETWEEN  SCANS  AND  THIS  SUBROUTINE 

6530  !  SHIFTS  THE  DATA  SO  THAT  THE  FIRST  DIODE  DATA  POINT  IS  MOVED 

6540  !  TO  THE  VERY  BEGINNING  OF  ITS  1124  BLOCK  IN  THE  OVERALL  ARRAY 

6550  !  .THE  FIRST  SET  IS  RIGHT  ON,  THE  NEXT  IS  ONE  OFF, THE  THIRD  IS 

6568  !  TWO  OFF  ,,SG  FORTH.  THIS  MAY  NOT  MATTER  WITH  OUR  RESOLUTION 

6570  !  AND  IS  PROBABLY  DUE  TO  THE  MEMORY  CARD  CYCLING  AT  THE  END  OF 

6580  !  SOME  SCANS.  SEE  THE  CIRCUIT  TIMING  DIAGRAM  IN  THE  THESIS. 

6590  SELECT  B 

6600  CASE  4  ! EXHAUST  DATA  HAS  4  SCANS 

6610  H=3 

6620  CASE  8  ! MOTOR  DATA  HAS  8  SCANS 


6630 

6640 

6650 

6660 

6670 

6680 

6690 

6700 

6710 

6720 

6731 

6741 

6750 

6760 

6770 


M=7 

END  SELECT 
FOR  X=1  TO  2 

IF  X=1  THEN  MAT  E=  Y1  I0PPERATES  ON  NO-PARTICLE  AND  PARTICLE  SETS 

IF  X=2  THEN  MAT  E-  Y2 
FOR  J=0  TO  M 

IF  H=7  AND  J<=3  THEN  6800  !One  4096  Block  Doesn't  Need  Shiftinq 
FOR  I=(J)*1024+1  TO  ( J+l )*1 024  ! BLOCKS  OF  1024 

IF  M=3  THEN  L=I+J  ! Array  B  Has  the  worst 

IPrsblen  With  Shifting  Data 

IF  M=7  THEN  L=I+1  ! Array  D  is  Allways  off  by  one 

IF'  ,J+1)»1024  THEN  L=(J+1)*1 024 
! JUST  TO  AVUID  PROGRAMMING  ERROR  AT  THE  END  OF  THE  ARRAY 
E(I  )=E(L )  ! THIS  SHIFTS  THE  DATA 

! DEPENDING  ON  'L', ARRIVED  AT 


1 1 


82 


V-V-Y- 
V --Y- 


>“  ■’ .  •*-*•*.  **«*"•  •*,  •*  .  •*.  «*,  •*  *  **,  •’ ,  /  *  ,  «*, 


!FY  OBSERVING  RAW  BATA 


6780 

6790  NEXT  I 

6800  NEXT  J 

6810  IF  K=1  THEN  HAT  Yl=  E 

6820  IF  K=2  THEN  HAT  T2=  E 

6830  NEXT  X 
6840  SUBEND 

6850  SUB  Store  !  THE  LARGE  ARRAY  CONTAINING  HULTIPLE  SCANS  HAS  BEEN  REDUCED 

6860  !  TO  A  HEAN  SCATTERING  PROFILE  BY  AVERAGING  AND  FILTERING.  IF 

6870  !  YOU  FEEL  THAT  STORING  THIS  DATA  IS  NECESSARY  THIS  ROUTINE 

6880  !  DOES  IT,  IT  SAVES  DISK  SPACE  TO  STORE  THE  RESULTS  THEN 

6890  !  ELIMINATE  THE  RAU  DATA  IF  YOU  FEEL  CONFIDENT  THAT  THE 

6900  !  REDUCTION  IS  THE  BEST  THAT  CAN  BE.  IN  OTHER  WORDS, 

6910  !  DO  NOT  PURGE  A  RAN  DATA  FILE  UNLESS  YOU  ARE  ABSOLUTELY  SURE 

6920  !  YOU  WON'T  WANT  TO  REDUCE  IT  AGAIN. 

6930  COH  /Hrdgaus/  Av2(»)  !Au2(*)  IS  THE  REDUCED  DATA 

6940  DIN  Dated:  1024)  BUFFER  ! BUFFER  USED  TO  TRANSFER  TO  DISK 

6950  HAT  Data=  Av2 
6960  PRINT  CHRSI12) 

6970  PRINT  USING  7 /////" 

6980  PRINT  *  A  SUGGESTED  METHOD  OF  NAMING  REDUCED  DATA  FILES  IS  AS  FOLLOWS’ 

6990  PRINT  ’  H - MOTOR  BEAM’ 

7000  PRINT  ’  X - EXHAUST  BEAM’ 

7010  PRINT  ’  C - CALIBRATION  ,  IF  NO  'C'  THEN  AN  ACTUAL  FIRING  IS  ASSUMED’ 

7020  PRINT  ’  G - GLASS  BEADS  ’ 

7030  PRINT  ’  IF  NO  'C'  THEN  'G'  STANDS  FOR  'GAP'  PROPELLANT’ 

7040  PRINT  ’  A - ALUMINUM  OXIDE’ 

7050  PRINT  ’  PPP - CHAMBER  PRESSURE  FOR  RUN  OR  OTHER. .. .(NOZZLE  TYPE)’ 

7060  PRINT  ’  PARTICLE  SIZE  FOR  CALIBRATION’ 

7070  PRINT  ’  MMDD - MONTH, DAY  OF  RUN  OR  CALIBRATION’ 

7080  PRINT  ’’ 

7090  PRINT  ’  EXAMPLE:  HCA125JN12’ 

7100  PRINT  ’’ 

7110  PRINT  ’  MOTOR  BEAM  CALIBRATION  USING  ALUMINUM  OXIDE  FROM  I  TO  25  MICRONS’ 
7120  PRINT  ’  ON  JUNE  12* 

7130  PRINT  ’  EXAMPLE:  XG550AU10’ 

7140  PRINT  ’’ 

7150  PRINT  ’  EXHAUST  BEAM  DATA  FOR  GAP  PROPELLANT  AT  550  psi  ON  AUGUST  10’ 

7160  PRINT  ’  PLACE  A  DISK  IN  THE  RIGHT  HAND  DRIVE’ 

7170  PRINT  ’  ENTER  THE  FILENAME  YOU  WISH  TO  USE  FOR  THIS  REDUCED  DATA’ 

7180  INPUT  At 

7190  CREATE  BOAT  At, 512, 16 

7200  ASSIGN  BDisk  TO  At 

7210  ASSIGN  BBuff  TO  BUFFER  Dataf*) 

7220  CONTROL  BBuff ,3jl ,8192,1  ! BUFFER  IS  FULL 

7230  TRANSFER  BBuff  TO  BDisk; COUNT  8192 
7240  WAIT  FOR  EOT  BDisk 
7250  ASSIGN  BBuff  TO  « 

7260  ASSIGN  BDisk  TO  * 

7270  SUBEND 


7280  SUB  Review(Avl(*)) 

7290  DI«  Datad:  1024)  BUFFER 

7300  PRINT  *  EACH  REDUCED  FILE  CONTAINS  ONE  SCAN  ONLY.  YOU  HUST  KNOW  IF  IT' 
7310  PRINT  ‘  IS  EXHAUST  ,  MOTOR  CAVITY  ,  OR  CALIBRATION  DATA.* 

7320  PRINT  ” 

7330  PRINT  1  THE  DISK  WITH  THE  REDUCED  FILE  SHOULD  BE  IN  THE  RIGHT-HAND  DRIVE1 
7340  PRINT  ■* 

7350  PRINT  1  ENTER  THE  FILENAME  OF  THE  REDUCED  DATA1 

7360  PRINT  *  TO  BE  REVIEWED1 

7371  PRINT  11 
7380  INPUT  A$ 

7390  ASSIGN  gDisk  TO  A$ 

7400  ASSIGN  BBoff  TO  BUFFER  Data!*) 

7410  CONTROL  gDisk, 5; 1 

7420  CONTROL  8Btiff,3;l, 0,1  (THIS  IS  AN  EMPTY  BUFFER 

7430  TRANSFER  gDisk  TO  SBuff (COUNT  8192  ! 1 024*8=NUHBER  OF  BYTES 

7440  WAIT  FOR  EOT  gDisk 
7450  ASSIGN  gDisk  TO  * 

7460  ASSIGN  gBuff  TO  » 

7470  HAT  Avl=  Data 
7480  SUBEND 

7490  SUB  Filter(A(*),Fil) 


7500 

DIH  E<1 : 1624) 

-v 

7510 

PRINT  CHRS112) 

-  li 

7520 

PRINT  USING  *////////* 

7530 

PRINT  1 

FILTERING  THE  SCATTERING  PROFILE1 

7540 

FOR  J=1  TO  Fil 

7550 

FOR  1=1  TO  1014 

7560 

B=I+1 

ITHIS  IS  A  Synetric  Hovina  Averaqe 

•*>V 

.  *  *  *  <( 

7570 

C=I+2 

•TYPE  OF  DIGITAL  FILTER .  EACH 

7580 

D=I+3 

IDATA  POINT  IS  EQUALLY  WEIGHTED 

— 1 

7590 

F=I+4 

f IN  THIS  CASE  BUT  THIS  CAN  BE 

7600 

G=I+5 

! CHANGED  IF  ONE  DETERMINES  THAT 

7610 

H=I+6 

'FEWER  POINTS  WITH  UNEQUAL  WEIGHTS 

7620 

0=1*7 

IWOULD  BE  FASTER  OR  GIVE  BETTER 

,*  ’/j 

7630 

L=I*8 

IRESULTS.  THIS  TYPE  OF  FILTER  WAS 

►  ‘  *“ 

7640 

N=I+9 

IUSED  SINCE  IT  INTRODUCES  NO  PHASE 

_  _  ^ 

7650 

P=I*10 

'LAG  (Anqtlar  Resolition) . 

7660 

E(G)=(A(I)+A(B)+A(C)+A(D)+A(F)+A!G)+A(H)+A(K)+A(L)+A(N)+A(P))/11 

tv 

7670 

NEXT  I 

7680 

HAT  A=  E 

\\VJ 

7690 

NEXT  J 

:  - 

7700 

SUBEND 

-  .  ^ 
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